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We investigate the universal orientation-dependent interactions between non-spherical colloidal 
particles immersed in a critical solvent by studying the instructive paradigm of a needle embedded 
in bounded two-dimensional Ising models at bulk criticality. For a needle in an Ising strip the 
interaction on mesoscopic scales depends on the width of the strip and the length, position, and 
orientation of the needle. By lattice Monte Carlo simulations we evaluate the free energy differ- 
ence between needle configurations being parallel and perpendicular to the strip. We concentrate 
on small but nonetheless mesoscopic needle lengths for which analytic predictions are available for 
comparison. All combinations of boundary conditions for the needles and boundaries are considered 
which belong to either the "normal" or the "ordinary" surface universality class, i.e., which induce 
local order or disorder, respectively. We also derive exact results for needles of arbitrary mesoscopic 
length, in particular for needles embedded in a half plane and oriented perpendicular to the cor- 
responding boundary as well as for needles embedded at the center line of a symmetric strip with 
parallel orientation. 



I. INTRODUCTION 

A solvent near its bulk critical point induces long- 
ranged forces between immersed mesoscopic particles 
which are "universal", i.e., independent of most micro- 
scopic details of the system [l-Q] . Apart from a few bulk 
properties of the solvent and from its local interactions 
with the surfaces of the immersed particles, these forces 
depend only on the geometry of the confinement of the 
critical fluctuations of the solvent imposed by the parti- 
cle surfaces, that is on the sizes, shapes, positions, and 
orientations of the particles. Such forces have been ob- 
served experimentally in film geometry for ^He [5|, |6[ 
and '^He/'*He mixtures [31 near the superfluid transition 
as well as for classical binary liquid mixtures near their 
demixing transition [8| and, directly, between a single 
spherical particle immersed in a binary liquid mixture 
near the critical demixing point and its confining pla- 
nar wall [3, [l3l- Due to the similarity of these forces 
with the Casimir forces arising in quantum electrody- 
namics [111, ll2| they are called critical Casimir forces. 

Nonspherical particles near other confinements, such 
as a planar wall, experience orientation dependent forces 
giving rise to critical Casimir torques [ij, [lj| . 

In the following we concentrate on critical solvents be- 
longing to the Ising universality class such as demixing 
classical binary liquid mixtures. In these systems the 
particle surfaces generically prefer one of the two com- 
ponents of the mixture, i.e., in Ising language one of the 
two directions (+ or — ) of the order parameter is pre- 
ferred. This preference is captured in terms of surface 
fields the strength of which under renormalization group 
fiow attains ±cx), corresponding to fixed surface spins 
and denoted by ± boundary conditions. However, using 



suitable surface preparation one can suppress this pref- 
erence, being left with a weakened tendency to demix 
near the particle surface [l5[ . These two types of surface 
universality classes [ig] are called "normal" (+/— ) and 
"ordinary" (O), respectively. 

Simple universal behavior arises in the scaling region 
where the ranges of the interactions between the order- 
ing degrees of freedom other than short-ranged ones and 
the lengths characterizing corrections to scaling and the 
crossover from less stable bulk and surface universality 
classes, respectively, are much smaller than the increas- 
ing correlation length of the bulk solvent upon approach- 
ing criticality, than the distances between the particles, 
and than the lengths characterizing the sizes and shapes 
of the particles. In this region of a clear separation of 
"microscopic" and "mesoscopic" lengths, the forces and 
torques are given by universal scaling functions which, 
apart from the surface universality classes, only depend 
on ratios of the characteristic mesoscopic lengths. 

Here, our main interest is in the orientation depen- 
dent interaction of nonspherical particles with bound- 
aries of critical systems. For the interaction between 
a prolate uniaxial ellipsoid and a planar wall with H — h 
or H — boundary conditions at the corresponding sur- 
faces the complete universal scaling functions for the crit- 
ical Casimir forces, both for the disordered and ordered 
bulk phases, have been obtained within mean field the- 
ory |14l ] . Beyond mean field theory such scaling functions 
have been determined to a lesser degree of completeness. 
For nearby particles, i.e., for closest surface to surface 
distances much smaller than their size and radii of cur- 
vature, the interaction can be expressed in terms of the 
Casimir force in film geometry by using the Derjaguin 
approximation. The scaling function of the latter is avail- 



able for ++ and H — surfaces in three |T7HT9{ and two [2^ 
spatial dimensions 21| • For a small spherical [22[ or non- 
spherical [13, I23 particle, with a size much smaller than 
the bulk correlation length and the distances to the sur- 
faces of other particles or to boundaries, the interactions 
can be obtained from the so-called "small particle expan- 
sion" [2JI . This is one of the field theoretic operator ex- 
pansions for small objects the most well known of which 
is the product of two nearby operators first considered 
around 1970 by Wilson, Polyakov, and Kadanoff [H,!!!]. 
These expansions are to a certain extent reminiscent of 
the multipole expansion in electrostatics. 

In view of the twofold asymptotic condition, i.e., the 
particle being small on mesoscopic scales and large on 
the microscopic scale [2^, it is a nontrivial issue under 
which circumstances the results of the "small particle 
expansion" can be observed for a given actual system. 
As a first step to address this issue, here we investigate 
whether the results of the expansion can be observed in 
the two-dimensional Ising model on a square lattice with 
(ferromagnetic) couplings J^ > between nearest neigh- 
bor spins only and right at the bulk critical point. 

In our Monte Carlo simulations we consider a lattice of 
finite size with the form of a strip (or rectangle) compris- 
ing W rows and L columns with an embedded particle 
resembling a needle. The directions along the rows and 
columns define, respectively, the directions u and v par- 
allel and perpendicular to the strip of length [231 ^ ^^'^ 
width W. While in the u direction we impose periodic 
boundary conditions by means of couplings JJ between 
the first and last spin in each row, wc couple all of the 
spins in the lowest and uppermost row by means of in- 
teractions or ±i7 to spins which are fixed in the -I- di- 
rection and are located in outside neighboring rows (see 
Fig. [T|) . Thus in the v direction, the strip is bounded by 
free surfaces or surfaces to which a magnetic field ex ±J' 
is applied, i.e., by surfaces of "ordinary" (O) or "normal" 
(-1-/—) character. 

In order to generate an embedded needle with O or 
+/— boundary conditions and orientation parallel to the 
strip we remove couplings (break bonds) between two 
neighboring lattice rows or we fix spins within a single 
lattice row, respectively (see Fig. [T]). In the two cases 
we define the length D of the needle as the number [23] 
of broken bonds and of fixed spins, respectively, which 
we choose to be even and odd, respectively, in order that 
the needle centers x coincide with the center of an ele- 
mentary square and a vertex of the lattice, respectively. 
This allows us to "turn" the needles abruptly about their 
center by 90 degrees upon rearranging the broken bonds 
and fixed spins, leading to broken bonds between neigh- 
boring columns and fixed spins within a single column of 
the lattice, respectively. In order to be able to position 
the centers of O and +/— needles right at the midline of 
the strip (as shown in Fig. [T]) we choose W to be even 
and odd, respectively. 

In the simulation we calculate the free energy expense 
AF required to "turn" the needle about its center from 
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FIG. 1: Needle, with their centers denoted by x, embed- 
ded in a strip of a ferromagnetic Ising model on a square 
lattice comprising W rows and L columns [23 of fluctuating 
spins (empty circles) and periodic boundary condition in the 
u direction parallel to the rows. The two additional rows of 
spins fixed in the + direction (full circles) allow one to in- 
duce positive or negative magnetic fields at the surfaces of 
the strip by coupling them in a ferro- or antiferromagnetic 
way to the bottom and top rows of fiuctuating spins (see the 
main text), (a) and (b): Needle of D = 4 broken bonds in 
a strip with the number W of rows and the number L of 
columns equal to W = L — 8. (a) shows the case {0[0]0) of 
a strip with two free surfaces (no coupling to the fixed spins) 
while (b) shows the case (-|-[0]0) in which a ferromagnetic 
coupling to the fixed spins leads to a positive magnetic field 
at the lower surface, (c) and (d) : Needle of D = 5 spins fixed 
in the -|- direction (full circles) in a strip with W = 7 and 
L = 8. (c) shows the case (0[-|-]0) while (d) shows the case 
(-f [-|-]— ) with positive and negative magnetic fields induced at 
the lower and upper surfaces by ferro- and antiferromagnetic 
couplings, respectively. All the needles shown are oriented in 
the u direction and have their center x at the midline w = of 
the strip, half way between the strip surfaces, i.e., vm ~ for 
the vertical position of the center of the needle. This requires 
to choose W even (odd) for needles of broken bonds (of fixed 
spins) so that spins reside at half odd integer (integer) values 
of II. Needles oriented in the v direction are shown in Figs. |6] 
and [3 below. 



an alignment perpendicular to the strip to the parallel 
alignment. This is a measure for the effective torque act- 
ing on the needle. In line with the introductory remarks 
the free energy cost AF depends on the distances of the 
needle center from the strip surfaces. 

We put the origin of the (u, v) coordinate system at the 
midline of the lattice, i.e., for O and +/— needles at the 
center of an elementary square and at a vertex so that 
the lattice spins are located at half odd integer and in- 
teger values, respectively, of the coordinates (see Fig.[T|). 
With this choice the mirror symmetry of AF in a strip 



with equal boundary conditions i = j is described simply 
by its invariance if the coordinate v = vj^ of the needle 
center changes sign. Here i and j denote the boundary 
conditions at the surfaces at w = —W/2 and v = W/2, 
respectively. Figure [I] shows the special case in which the 
centers x of O and +/— needles coincide with the ori- 
gin, i.e., vn vanishes, while in the general case vjy takes 
integer values for both types of needles. 

It is useful to adopt a notation {i[h]j) which character- 
izes the boundary conditions i and j of the strip surfaces 
as well as the boundary condition h of the embedded nee- 
dle. For example, (+[+]—) denotes the case of a needle 
of spins fixed in the -I- direction which is embedded in a 
strip with outside couplings J^ and — ^J at the lower and 
upper surfaces inducing strip surfaces i ~ + and j = —, 
respectively, as shown in Fig. [TJd). 

For completeness we consider also the case of a strip 
with periodic boundary conditions in both u and v di- 
rections by coupling, in both the W rows and in the L 
columns, the corresponding end spins to each other with 
strength J^. This is a strip without surfaces and is equiv- 
alent to a square lattice on a torus. 

If D, W, L, and the closest distance between the nee- 
dle and the strip boundaries (corresponding to gn and 
a< in Figs.^a) and[21[b), respectively) are "sufficiently" 
large on the scale of the lattice constant, the free en- 
ergy cost AF in the lattice model is expected to dis- 
play universal scaling behavior and, as mentioned above, 
AF/{kBT) will depend only on the universality classes 
{i[h]j) of the boundaries and the needle and on three in- 
dependent ratios of D, W, L, and vj\[ . In this case one 
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FIG. 2: Continuum [231 description of the geometry of a nee- 
dle in a strip. In the strip of length L and width W the 
embedded needle of length D is oriented parallel (a) and per- 
pendicular (b) to the strip boundaries ai v — iizW/2. The 
needle center x is located ai u — un ~ and v — vn < 
at a distance un = vn + W/2 > from the lower boundary. 
In the perpendicular orientation the closer end of the needle 
has a distance a< = ajv — D/2 > (and the farther end a 
distance a> = a< -I- D) from the lower boundary. 



can adopt a mesoscopic continuum description with sharp 
strip boundaries at w = ztW/2 since lengths only differ- 
ing by approximately a lattice constant can be identified. 
Likewise, different microscopic definitions of the length 
of the needles, such as the number Z3 — 1 of bonds be- 
tween the fixed spins of a -f needle (instead of its D fixed 
spins), lead to the same mesoscopic length D. Figure [2] 
shows the various lengths characterizing the geometry of 
a needle in a strip within the continuum description. 

For small mesoscopic needles and a long strip, i.e., for 
-D <C W^ <C i, the universal behavior can be predicted 
by means of the "small particle expansion" and, as men- 
tioned earlier, it is one of our main goals to quantitatively 
investigate whether present Monte Carlo simulations can 
access this regime. 

For a needle perpendicular to the boundary of a half 
plane or embedded at the midline of a symmetric (i = j) 
strip of infinite length, we derive the analytic form of the 
universal scaling behavior of the critical Casimir forces 
in the complete range of mesoscopic needle lengths D, 
increasing from small to large. This aspect is interesting 
in its own right and also allows us to estimate the range 
of validity of the "small needle expansion" . 

The predictions of the "small particle expansion" for 
AF in critical Ising strips and our analytic results for nee- 
dles of arbitrary length are presented in Sees. lU and IIIIl 
respectively. In Sec. II VI we describe the lattice model in 
more detail and explain how one obtains via Monte Carlo 
simulations the results for AF which arc compared with 
the analytic scaling predictions in Sec. |V] . Our results 
are summarized in Sec. |VIl In Appendices E] and |B] we 
present the input material and the derivations which are 
necessary in order to obtain the predictions and results of 
Sees, im and mil respectively. For the convenience of the 
reader we present a glossary of our symbols and notations 
in Appendix C. 

Besides for critical systems [ij, [Ij] the orientation de- 
pendent interaction between a wall and a nonsphcrical 
mesoscopic object (such as an ellipsoid, a semi-infinite 
plate, or a sphero cylinder ) ha s been studied also for quan- 
tum electrodynamics [2^ [23 and for purely entropic sys- 
tems in which it is induced by spherical so-called deple- 
tion agents with hard body interaction only [SOf or by free 
nonadsorbing polymer chains [23|. As for correlation- 
induced forces we mention also the attractive effective 
force generated by needles (rigid rods) acting as deple- 
tion agents [3l| and the repulsive force generated by 
a nonadsorbing polymer chain grafted to the tip of a 
(model-) atomic force microscope near a repulsive wall 
j32| . Finally we note that the critical Casimir force be- 
tween inclusions in the two-dimensional Ising model has 
been suggested as a possible mechanism for the presence 
of long-ranged forces between membrane bound proteins 
[33|, which arc typically non-circular. 



II. SMALL NEEDLE EXPANSION 

Here we specialize the "small particle expansion" for 
non-spherical particles (see Ref. [ij]) to the present case 
of a needle embedded in the two-dimensional Ising model 
and apply it to the geometry of a needle in a strip. 

We consider a needle of small mesoscopic length D, 
centered at rjv, and pointing along the unit vector n. 
Inserting the needle into the d = 2 Ising model at its 
critical point [34| changes the Boltzmann weight of the 
corresponding field theory by a factor exp(— (JH) which 
can be represented by the operator [2^ expansion [Tj, [23 
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are the isotropic and anisotropic contributions, respec- 
tively. Here O = (j) is the order-parameter-density op- 
erator and O = e is the difference of the energy-density 
[3a | operator e and its average (e)buik in the unbounded 
plane (bulk) at bulk criticality (so that their bulk aver- 
ages ('/')|-mii; and (e)buik vanish at the bulk critical point). 
They arc normalized such that the bulk two-point corre- 
lation functions [341 are ISol 
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where (in d = 2) x^ = 1/8 and x,, ~ 1 are their scaling di- 
mensions [37| . The affiliation of the surface of the needle 
to the "ordinary" {h ~ O) or to the "normal" {h = + or 
/i = — ) surface universality class enters into Eqs. (|2.ip - 
via the universal half plane amplitude A^^ . The 

of O in the half 



latter is the amplitude of the profile j34 

plane, (C')haif piano = -4^0"^°, as function of the meso- 
scopic distance a from the boundary line of type h and 
is given by [33 
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The properties A\ = and A\ 
plitudes of the order parameter profile are obvious from 
the up-down symmetry of the Ising spins [3J] . The en- 
ergy density profile increases (decreases) upon approach- 
ing an "ordinary" ( "normal" ) boundary where the Ising 
spins are more disordered (more ordered) than in the 
bulk [35| . This implies the positive (negative) sign of the 
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to duality propertie s l38l . |39[ . In Eq. (|2.3p the contribu- 
tion from the stress [331 tensor Tki is the same [401 f^^' ^1^ 
needle types /i = O, +, — . 

The ellipses in Eqs. (|2.2[) and (|2J}[) stand for contribu- 
tions of fourth order and higher [4l| in the small meso- 
scopic length D. The common factor of proportionality 
on the rhs of Eq. (|2.ip is given [33 by (exp(— (5'H))buik 
because the bulk averages of Sj and Sa vanish at the 
bulk critical point. 

For a two-dimensional Ising model with boundaries the 
insertion free energy of the needle depends on its posi- 
tion and orientation with respect to these boundaries. 
Removing the needle from the bulk model and inserting 
it in the bounded one at a distance from the bound- 
aries much larger than D changes the free energy by 
F{rN,n) = -ksT \n{{l + Si + Sa))bm where BM denotes 
the bounded model in the absence of the needle. We shall 
concentrate on the geometry of a needle in a strip as de- 
scribed in the Introduction (see Fig. [5]) for which BM = 
ST is the strip with boundaries {i,j) in the absence of 
the needle and correspondingly use the notation [23| 
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for the center and direction vectors of the needle with the 
components parallel and perpendicular to the strip. The 
density averages (C')st at tn which enter {Sj + Sa)st 
are, within this model, independent of un and in the 
scaling region given by 
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where / are universal scaling functions and ./q'' follows 
from /q'"' upon replacing vpj/W by —vn/W. 

Due to the continuity equation of the stress tensor [33 
its average {Tki{rN))g^' in the strip is independent of 
both UN and vn and follows from the universal, scale- 
free, and shape-dependent contribution [42[ ^^.^'{L/W) 
to the free energy Fg^ per ksT of the strip ST without 
the needle: 
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Like in Eq. (|2.6p . here || and _L denote directions parallel 
and perpendicular to the u axis of the strip, 

X,j{W/L) = id/dS)^l!T!\S), 6 = L/W, (2.9) 



and the off-diagonal components (rL||)sT = (T'||_l)st 
vanish by symmetry. Cardy [43| has obtained an ex- 

pheit form for all functions ^g!^"* (iS). While an ex- 
tensive discussion of Aij(W/L) is deferred to Ap- 
pendix IA21 here we mention a few basic properties. 
J .J are symmetric in i and j. 
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for (i,j) = [(0,0),(+,+),(+,-),(+,0)]. While due to 
the (-1- -H' — ) symmetry the equalities A+_|.(1/(S) = 
A_,_(l/(5) and A+,o(l/(5) = A_.o(l/(5) hold for arbi- 
trary 1/(5, the equality Aqo = ^+,+ holds only for in- 
finitely long strips [4J], i.e., 1/6 = 0. 

For the free energies i^n and F± associated with remov- 
ing the needle from the bulk system and inserting it in 
the strip with its center at r^r and with its orientation 
n parallel and perpendicular, respectively, to the u axis, 
Eqs. dm])- (1131) and dill) yield 
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follow from Si and Sa in Eqs. (|2.2p and 
tively. This implies the expression 

AF = -kBTln{ZJZ±) 



...(2.14) 
respec- 

(2.15) 



for the free energy required to turn the needle about its 
center from its orientation perpendicular to the u axis of 
the strip to the parallel orientation. 

In the expansions in Eqs. (|2.13p and p.l4p a term 
(X D"^ is of the order of {D/W)'^ near the strip cen- 
ter, where \vn\ "C W, and of the order of (D/qn)'^ near 
the strip boundaries, where a^ = {W/2) — \vn\ ^ W. 



These expansions for the partition functions Z\\ and Z}_ 
are useful if D/W and D/on are much smaller than 1. 
However, expanding the free energies Fn, F±, and AF 
in powers of D is not always useful for the compari- 
son with our simulation data. While D/W <^ 1 and 
D/qn <C 1 can readily be achieved in the MC simula- 
tions for a needle with a mesoscopic length D of many 
lattice constants, at present it seems to be unrealistic to 
achieve {D/W^^^ < 1 and {D/onY/^ < 1. Thus we 
expand the logarithm in Eq. (j2.15p in terms of powers of 
D only in those cases in which the power Z)^* = D^^^ 
does not appear [43 in Eq. (|2.13p . These are the cases 
with an "ordinary" needle h = O, for which Aj^ van- 
ishes, and with a "normal" needle in a strip with two 
"ordinary" boundaries i ~ j = O, for which the den- 
sity profile in Eq. (|2.7p vanishes for the order parameter 
O = 0. Expanding the logarithm in these cases yields 

AF^AFi + AFr.i + ... (2.16) 

with the leading contribution 



and the next-to-leading contribution 
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The ellipses in Eq. (|2.16p are of order _D*. 

For completeness we also consider a double periodic 
rectangle or strip with periodic boundary conditions in 
both the u and w-directions, i.e., the surface of a torus. In 
this boundary-free case at Tc the average of vanishes in 
the strip, that of e is independent of both u and v, and 
for the average of the stress tensor Eq. (|2.8I) again ap- 
plies with Ai,j{W/L) replaced by a function Ap{W/ L) = 

{d/d5)(^fr^ {5) with the infinite strip hmit Ap(0) = Ap = 
— 7r/12 (see Subsec. FV Al and Appendix I A II for more de- 
tails). For the double periodic strip Eqs. (PT^ - (PJ5)) 
also apply if A,j(VF/L) and f^''^\vN /W,W/L) are re- 
placed by the corresponding stress amplitude Ap(W^/L) 
and the energy density scaling function /i ' (W/L) which 
is independent oi vn/W. 

Explicit expressions for the functions A{W/L) are 
known |43j for all types of strips considered here. Con- 
cerning the scaling functions fo the dependence on the 
aspect ratio W/L is known [4y, |43] for /^ while in the 

presence of boundaries the functions .f^n ''' (vn /W, W/L) 
are, to the best of our knowledge, known only |48| for 
W/L — 0, i.e., for strips of infinite [43 length L = 00. 
For the convenience of the reader we collect these expres- 
sions in Appendix [Xl 



For the special case in which the distance ajv = Vn + 
W/2 of the center of the small needle of type h from the 
boundary of type i at w = —W/2 (see Fig. [2]) is much 
smaller than the width of the infinitely long strip, i.e., 
for the case D <C ajv and W/a^, L/un — >■ oo, the above 
expressions for the needle in the strip with free energies 
F = i^(*[''lj) reduce to those for the needle in the half 
plane v + W/2 = a > with free energies F = i^(*[''l) 
where Eqs. (12.13^ and (|2.14p are replaced by 



1,7 — 2_^ ■^o ■^o 



0=4>,e 



^xc 



1+^O.q2+. 
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and Eqs. (PTf)) and ((^1^ by 

^^' 0, $.^=^w^w?^3 



knT 



knT 



where 



D 

2aN 



(2.19) 



(2.20) 



(2.21) 



As mentioned above these relations only apply if < ?? <C 
1. 



III. NEEDLES OF ARBITRARY LENGTH 

The "small needle expansion" is valid if the mesoscopic 
length D of the needle is "small" compared to the other 
mesoscopic lengths of the system, i.e., within the present 
model much smaller than the width and length of the 
strip and the distance of the needle from the two bound- 
aries. Given the limited set of operators appearing in 
Eqs. (|2.2p and (j2.3|) . it is an open issue what in this con- 
text "small" means quantitatively. In order to get a clue, 
in this section we study a few situations in which results 
for an arbitrary mesoscopic needle length D can be ob- 
tained. These results, which we derive in Appendix [B] 
are interesting also in their own right. 

(i) Needle in half-plane 

We consider a needle of universality class h embedded 
in the half plane, oriented perpendicular to the boundary 
line of surface universality class i, and with its center 
located at a distance oat from the boundary. Here we 
also introduce the distance a< = ajy — {D/2) of the closer 
end of the needle from the boundary (compare Fig. [^b)) 
so that one has 



d. 



1 



D 



2aN l + {2a</D) 



(3.1) 



for the length ratio 'd defined in Eq. ((2?2T|) . Note 
that d tends to and 1 for a small (or distant) nee- 
dle with D -^ a^ and a long (or close) needle with 
a< ^ D, respectively. Both D and a< are assumed to 



be mesoscopic lengths, i.e., both are large on the mi- 
croscopic scale. We are interested in the free energy 
Fl = F^/ ~ ^BTf^f^M) required to insert the needle 
from the bulk into the half plane, and in the correspond- 



ing Casimir force —{d/dapf)F\_ 



{t[h]} 



<d/da^)F 



im) 



[ksT / aN)g-/7f^Ad) with [50] universal scaling functions 
/j[^i and gi\h]- The above partial derivatives are taken at 
fixed needle length D. The force pushes the needle away 
from (towards) the boundary iigi[h] is positive (negative). 



Due to symmetries the identities F 



(+[+]) 



F 



(-[-]) 



F^+m ^ pi^-m ^^^ pirn) 



fI'^'*'^ hold so that ex- 



changing the surface universality classes of the needle and 
of the boundary leaves the free energy unchanged. These 
follow from the (+ ^-^^ — ) and the duality [51| symmetries 
of the Ising model and are consistent with the symme- 
tries of the corresponding small needle expression which 
follows from Eqs. (piT|) . (j^l^ . and ([^T^ . 

(ia) For i ~ h = O the effective interaction has the 
form (see the paragraph containing Eqs. (|B2p - (|B4p ') 



l^{0[0]) 

_± 

knT 



foio]W 



-2-Mn 



(1 + 1^)^ 

(1 - ^)3 



(3.2) 



and is attractive within the entire range < i? < 1. 
For z? ^ 1, Eq. (|3.2p is in agreement with the cor- 
responding resuh FJ^'^'VC^s^) = -^/4 + ^^/'^^ - 
')9'^/12 + 0(1?^) of the small needle expansion which fol- 
lows from Eqs. ((2lT|) . ([2l^ . and (|2l^ . For ■&/■!, 
Ff^'^^'^/keT -)■ -2-5 X 31n[(22/3D)/a<] shows a loga- 
rithmic dependence. The logarithmic divergence of the 
free energy for D/a^ — > 00 is related to the long-ranged 
behavior ex l/a< of the Casimir force for a needle of infi- 
nite length D — 00 which is addressed in Eq. p.9p below; 
its integral diverges for a< — > 00 (while for a needle of 
finite length D, the force decays more rapidly than l/a< 
as a< increases beyond D so that the integral is finite), 
(ib) For i = -\-, h = O the interaction is always repul- 
sive and has the form 



F 



i+[o]) 



knT 



'In 



[l + 'df 



(3.3) 



i.e., the same form as Eq. p.2p but with ■& replaced by 
~'d (see Appendix IB ip . 

(ic) For the Casimir forces, in Appcndix lB 2l we find ex- 
pressions which apply to arbitrary combinations (i, h) of 
the universality classes. For later convenience we present 
them here in two forms 



d F 



(m) 



daN ksT 



A0[0]) 



uh{'^) 




(3.4) 



which are equivalent due to Eq. (13. 2|) and the definition 



The dependence of pih on i? is given by 



(3.5) 



witli the variable 1/(5 replaced by the function 

l/S = K*{d)/['iK{'d)] (3.7) 



of d. Here 



K*{§)=K{i}), i9 = Vl - ??2 , 



(3.8) 



and ii' is the complete elliptic integral function (see 
Eqs. (8.113.1) and (8.113.3) in Ref. [S^]). The quantity 
Ai,,i(l/(5) in Eq. ^^ is the well-studied ^ Casimir (or 
stress tensor) amplitude for a strip ST without needle, 
with boundaries («, h), and finite aspect ratio W/L ~ l/S 
which has been introduced in Eqs. (|2.8p and p.9p and is 
given explicitly in Appendix IA2I 

^^''''1^ addressed above Eq. ^^ 
(13.81) by the symmetries of 



The symmetries of F]^ 
are reflected in Eqs. (13.41 



Ai,j{l/S) discussed below Eq. (j^ . 

For the special cases (i, h) = {O, O) and (+, O) we 
have checked that Eqs. p.4p - p.8p are consistent with the 
simple expressions given in Eqs. p.2p and p.3p . In par- 
ticular, fo.o('?) vanishes for all -d. 

The present results for a needle embedded in a half 
plane involve, via Ai^ft(l/(5), knowledge about the strip 
ST without needle but with finite aspect ratio because 
the two geometries are related by a conformal mapping 
as explained in Appendix IB 21 

Varying •& from to 1, for all combinations {i,h) 
Eqs. p.4p - p.8p provide the complete crossover of the 
force ~{d/daN)Fj^ /ksT as the length D and the posi- 
tion (determined by the distances an or a<) of the needle 
in the half plane change from small and remote from the 
boundary to large and close to the boundary. Equation 
p.7p tells us that the corresponding change in the ge- 
ometry of the strip ST is from remote [5 \ 0) to close 
((5 /* oo) strip boundaries, as expected. 

In the small needle limit ^ \, the quantity fj jj van- 
ishes (compare Eq. (|A13|) ) and the expansion for small d 
is provided b y Eqs . (jM)) . (|A10|) . and (|A12|) and the form 
of a in Eq. (jB22p . We check in the paragraph contain- 
ing Eq. (jB26p that this expansion is consistent with the 
small needle expansion in Eqs. (|2.1ip . (|2.12l) . and (|2.19p . 

In the long needle limit d /^ 1 one has fiji ^• 
16[A,,i(0)-A,,,,(0)]/7r with Ai,/i(0) given below Eq.' ([2J)) 
so that [531 



-(a/9ajv)F|^"''VfcBT ^ [-3, 
(z[/.]) = [(0[0]), (+[+]),(- 



-3, 61, 5]/(32a<), 
+]),(+[0])]. (3.9) 



The behavior ex a< of the force per kBT acting on the 
needle of infinite length D = oo follows easily [53[ from 



comparing its scaling dimension with that of the only 
remaining mesoscopic variable a< . Unlike the expansion 
of the force for a small needle in which fractional powers 
such as ajf (Z?/ajv)^'* may occur, the expansion around 
the long needle limit involves only contributions for which 
Eq. p.9p is multiplied with positive integer powers of 
a^/D = a (compare Eqs. (|A7p and (jAlip . the relation 
between k and d in Eq. (|B22p . and the relation ?9^ = 
4a(l + a)/(l + 2a)2 which follows from Eq. ((3?T|) ). 
In the limit d /^ 1 the free energy per ksT tends to 



F^l^'^^/iksT) 



ft, 



[h\ 



^ ([-3, -3, 61, 5]/32) ln(C,, J(l - §)) , 
(z[/.]) = [(0[0]),(+H), (-[+]), (+[0])] 



(3.10) 



where 1/(1 — ^) — > I?/(2a<) in terms of a<. Here Ci^h 
are numbers, in particular Cq.o = 2^/-^ and C+,o = 
2^/5 (see Eqs. ^^ and ^^). This logarithmic be- 
havior of -FJ* /{^bT) for the long perpendicular nee- 
dle should be compared with the power law behavior 
F|| /(ksT) -^ AihD/a^ of the free energy for a long 
needle aligned at a small distance a< parallel to the 
boundary which is obtained by turning the perpendic- 
ular needle about that end which is closer to the bound- 
ary. For a particle of circular shape [22| with a radius 
R much larger than the distance a< between the clos- 
est points of the circle and the boundary, the interaction 
free energy F^'™/{kBT) -^ 7rA,,ft(2i?/a<)i/2 exhibits, 
as expected intuitively, a power law exponent the value of 
which lies in between those corresponding to perpendic- 
ular and parallel needles. Unlike the long perpendicular 
needle, for the long parallel needle and the large circle 
the exact results of the interaction free energies quoted 
above are found already within the Derjaguin approxi- 
mation and exhibit an amplitude proportional to Ai^h 
which depends on {i,h). For all three particles (the per- 
pendicular and parallel needles as well as the circle) the 
particle-boundary interaction is attractive (repulsive) for 
i = h {i ^ h). 

Figure [3j a) shows a comparison between the exact re- 

suh for F\^^^^^/{kBT) = ft[+]W which follows from 
Eq. (|3.4p upon integration [50| and its "small needle ap- 
proximation" ("sua") (compare Eqs. (|B27)) and (|B33|). 
The approximation reproduces the exact result very well 
for d <do with 0.3 <^o< 0.4 while it fails for d > 0.6 
where it shows an unphysical maximum. For the strip sit- 
uation (— [+]j) with the center of the perpendicular nee- 
dle at vn, the question arises down to which proximity 
of the -|- needle to the lower — boundary of the strip one 
can expect the "small needle approximation" to be rea- 
sonably valid. Identifying the distance ajv = F)/{2'd) in 
the half plane with the distance W/2 — \vn\ = W/2 + vn 
in the strip suggests the approximation for F 



{-[+]]) 

± 



to 



be valid if \vn\/W is smaller than [1 - DI{W^o)]/2. For 
0.3 < ^0 ^ 0.4 and the value D/W = 21/101 used in our 
MC simulations (see, c.f., Subsec. FV Cp this rough argu- 
ment predicts that IiiatI/M^ should be kept smaller than 
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FIG. 3: Free energy cost per fcsT, /Ai, to transfer a needle of length D and universality class h from the unbounded plane 
(bulk) to the half plane, with an orientation which is perpendicular to the boundary of the half plane of universality class 
i. The dependence on i? = D/{2aN), with ajv as the distance of the needle center from the boundary, is shown for various 
combinations {i[h\). The exact results for arbitrary D/2 < aN are compared with the "small needle approximation" ("sna", 
Eq. (|B33p ). which is valid for small i?, and with the long needle behavior for i9 — >■ 1 as given by Eq. (|3.10|l . The interaction 
between the boundary and the needle is attractive (repulsive) for equal (different) universality classes i = h {i ^ h) shown in 
(b) and (c) ((a) and (d)). 



about 1/6 or 1/4, i.e., vn should be kept rather close to 
the midline v = 0. For the anisotropy Ai^'^I+l-'^ of the 
free energy of the needle in the strip there is even more 
uncertainty concerning the validity of the "sna" because 
the corresponding comparison for Fu is lacking, i.e., 

the counterpart of Fig. ^a.) for a needle parallel to the 
boundary of the half space is missing. 

Figures[3l[b)-(d) show the corresponding results for the 
free energies fcsr/Ai with {i,h) = {0,0), (+,+), and 

(+, O). For i = h and i ^ h, f,j\-f^i is negative and positive, 
respectively. 

(ii) Needle in infinitely long symmetric strip 

Here we consider a needle of universality class h and 
with parallel orientation || at the center line u = of 
an infinitely long (i,i) strip of width W. We determine 



the free energy cost F\\ = F, 



i[h]i 



keTl 



l^m and 



the corresponding disjoining force — ((9/9VF)F|| '^' = 

{kBT/W)g[^^{6) between the two boundary lines of the 
strip upon inserting the needle from the bulk with par- 
allel orientation. Here, the partial derivative is taken 
with the needle length D kept fixed. The disjoining force 
increases (decreases) if giihU is positive (negative). The 

{i,h) symmetries of i^|| are those of Aiji as can be 

inferred from E g. (I3.15P below. The results depend on 
the length ratio [5J| 



e^ 



W 



(3.11) 



and provide the crossover from the small needle behavior 
for ^ 1, consistent with the "small needle expansion", 
to the long needle behavior for 6* ^ 1. In the latter case 



one has 



F, 



ilh]i) 



kuT 






-^ D 



A,,^(0) A,,,(0) 



{W/2) W 

= e[-l, -1, 31, 3]/16, 
(z, h) = [{O, O), (+, +), (-,+), (+, O)] (3.12) 

because upon inserting a long parallel needle with its 
center at m = w = transforms, within the interval 
-D/2 < u < D/2, the {i,i) strip of width W into two 
independent adjacent (i, h) strips of width W/2 each. 

(iia) For i = h = O the free energy to insert the needle 
from the bulk reads 



F, 



{0[0]0) 



keT 
(iib) For i = +, h 



1 sinh( 

— m 

16 9 



O wc find 



F, 



(+[o]+) 



knT 



1 , sinh6' 
16^"^ 



(3.13) 



(3.14) 



which is Eq. p.l3|) with 6 replaced by —9. Equa- 
tions (|3.13p and (|3.14p are derived in Appendix IB II 
(iic) For arbitrary {i,h) we find (see Appendix IB 2[) 



^(i[/i]i 



{010]0) 




1 9 

ml' 
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2Ayt^'}^lP--^'\ 



(3.15) 



where fi^h{t) and Pi^h{t) are the functions Ti^h{^) and 
Pi,h{'&) from Eqs. p.5|) and (|3.6p . respectively, evaluated 
at I? = i. Equations (|3.13p and (|3.14p as well as Eq. 
p.lSp are consistent with the "small needle expansion" 
for small 9 (see Eqs. ((B29| ') and with Eq. ([XT^ for large 
9. 

Figure |4] shows a comparison between the exact re- 



A-[+]~). 



which follows 



suh for F|p'^'~V(fcBr) = /V-]_(6') 
from Eq. p.lSp via integration [50| and its "small nee- 
dle approximation" (compare Eqs. (jO)) . (|B30p . (|B32p . 
and (|B34l) ). For the value 6* = 2l7r/ 101 ^ 0.653 as used 
in our MC simulations (see Subsec. IV Cp there is good 
agreement. However, we have no such comparison for 



F 



{i[h\i) 



■II 



All four possible independent free energies kBTfL-,^{9) 
for needle insertion are collected in Fig. [SJ They are 
negative (positive) for i ^ h {i ^ h). In this double 
logarithmic plot the power laws as obtained for small 9 
from the "small needle expansion" and for large 9 from 
the long needle expression (Eq. p.l2p ) appear as straight 
lines. 

(iii) Periodic strip 

Finally we consider an infinitely long strip of width W 
with periodic boundary condition in v direction. In this 




FIG. 4: Free energy cost per kBT, /_r , , , to remove a needle 
of universality class -I- and length D from the bulk and in- 
sert it along the midline of an infinitely long strip with (— , — ) 
boundary conditions and width W as a function of8 = ttD/W 
(see Ref. [50| and Eq. (|3.15|l '). The exact result for arbitrary 
length ratio 9 is compared with the "small needle approxi- 
mation" ("sna", Eq. (|B34p '). valid for small 9, and with the 
long needle behavior according to Eq. (|3.12p . Also shown are 
the first term, (26')^'*, and the sum of the first three terms 
[29^/^ + 6i^/*/2^/* + {29f'^/2,, denoted as "snaeS", in the 
expansion of /_r , i_ in terms of powers of 9 for 9^1 (see 
Eqs. (|B29P and (|B32p ). Expanding /ir , i_ proves to be less 

successful than the "sna" within which exp( — /_,,_) is ex- 
panded and which agrees with the exact result up to ^ 1. 
The first term (dash-double-dotted green fine) approaches the 
exact result (full red line) for ^ ^- only very slowly. 



case the free energy cost to insert a needle parallel to the 
strip reads (see Appendix IB ip 



fcflT 



1 sinh ( 

— m 

16 9 



(3.16) 



This result holds for any needle universality class h = 
O, +, — and it is valid for an arbitrary ratio 9 = ttD/W. 
Equation p.l6p should be compared with Eq. p.l3p . For 
61 > 1 the rhs of Eq. ([3T6l) tends to 6'/16 = Dx[Ah,h{0)- 
Ap{0)]/W. This is consistent with the periodic strip 
of width W being transformed, by inserting a parallel 
needle, to a (/i, h) strip of width W within a w-interval of 
length D. Consistent with the "small needle expansion" , 
for < 1 the rhs of Eq. ([XTC)) tends to 9"^/ 96 + (D{9^). 



IV. FREE ENERGY ANISOTROPY FROM 
MONTE CARLO SIMULATIONS 

The anisotropic shape of a bounded critical system 
induces orientation dependent properties for embedded 
non-spherical particles. In Sec. |TT] the operator expan- 
sion has provided the asymptotic scaling properties for a 
"small" but "mesoscopic" particle in a "large" system. In 
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FIG. 5: Same as Fig. U for all four independent combina- 
tions {i[h]i) of a needle of universality class h placed along 
the midline of an (i, i) strip as function oi 6 = irD/W. 
Full, dash-dotted, and dotted lines correspond to the ex- 
act results, the "sna", and Eq. (|3.12p . respectively. For 
h = i {h ^ i) /jLij is negative (positive), decreasing (in- 
creasing) with increasing 6, and implies a negative (posi- 
tive) contribution -awF|f''''''' = {kBT/W)edefl^]i to the 
disjoining force. Thus the attractive Casimir interaction 
(see Eq. (|2.8p ') with the negative universal disjoining force 
-dwksT^'^ST^'' = Ai,ikBTL/W^ = -{■K/A8)kBTL/W^ be- 
tween the (i, i) boundaries of the long strip without the nee- 
dle becomes even more attractive (less attractive) due to the 
presence of the needle with h — i (h ^ i). 



the present section we address the issue to which extent 
these asymptotic properties capture the actual behavior 
in specific critical model systems. Concerning the needles 
studied here, we want to check whether the asymptotic 
predictions of Sec. [II] for A_F can already be observed 
within a lattice model with numerous, but not too many 
rows and columns so that the model is amenable to sim- 
ulations. In this section we describe how to set up the 
corresponding Monte Carlo simulations and to calculate 
AF. In Sec. |V] we compare these simulation data with 
the corresponding analytic predictions. This allows us to 
judge both the achievements and the limitations of the 
"small particle approximation" . 



A. Model 

For the simulation we use the lattice version of the Ising 
strip described in the Introduction and shown in Fig. 
1 . The implementation of the double periodic boundary 
conditions is obvious. Beyond that, here we describe in 
more detail strips with boundaries {i,j). In this case the 
lattice Hamiltonian Hst for strips without a needle reads 
(see Sec. HJ: 



'Hst/J — — 2_^ Ju,u;u'. 



'^i 7 ,Su-{W+l)/2Su-{W-l)/2 
u 

~^J Z^^u,{W+l)/2Su,{W-l)/2 (4.1) 



with J7 > and Ju.v-.u' .v' = 1 for nearest neighbors (de- 
noted by {u,v;u' ,v')) and zero otherwise. The fluctu- 
ating Ising spins s = ±1 reside on the vertices of the 
square lattice consisting of W rows and L columns with 
periodic boundary conditions in the u direction. The 
two last terms generate "ordinary" or "normal" (+/— ) 
boundaries near v = ~W/2 and v = W/2 by means of 
fixed values Su,-{w+i)/2 — Sm,(w+i)/2 = 1 for the spins 
in the two additional outside rows and by choosing in- 



i(i) 



i(i) 



dependently AJ ' and A^- ' equal to ("ordinary") or ±1 
(-1-/— "normal"). This serves as a microscopic realization 
of all the pairs {i,j) of strip boundary types considered 
in Sec. HIl which exhaust all possible surface universality 
classes in d = 2 (up to the "extraordinary" boundary 
type |38l . |44| corresponding to infinitely strong couplings 
between surface spins). 

For reasons given in the Introduction and shown in 
Fig. 1, for "ordinary" ("normal") needles the number 
W of rows is taken to be even 2, 4, 6, ... (odd 1, 
3, 5, ...) and the components of the lattice vertices 
u, V are half odd integers ±1/2, ±3/2, ±5/2, ... (integers 
0, ±1, ±2, ...). The components u = un and v = vn of 
the needle centers x are integers for both types of nee- 
dles. 

In order to be able to compare the Monte Carlo data 
with the results of Sec. II we consider the system at 
its bulk critical point T^ with J/iksTc) = \ ln(v^ + 
1) where J is the coupling constant scaled out of "Hst 
(Eq. gH). 

Inserting a needle of class /i = O or ft, = ± into strips 
amounts to appropriately removing bonds or fixing spins 
in accordance with Fig. [T] via additional terms in the 
Hamiltonian. For given surface universality classes (i,j) 
of the strip, universality class ft, needle length D, and the 
coordinate vm of the center of the needle we introduce the 
notation 



n = Ho^ HsT 



y^ih) 



and 



H^Hi=HsT + H 



(h) 



(4.2) 



(4.3) 



^U^V '-*U' ^v' 



for the total lattice Hamiltonian H corresponding to the 
needle being oriented perpendicular and parallel, respec- 
tively, to the u axis. The explicit forms of "H^ and Ti,, 
will be given in Eqs. (g^, (gS]), and (|IlT|) - (|il^ below. 

B. Numerical algorithm 

1. Coupling parameter approach 

We consider two systems with the same configura- 
tional space C (i.e., number and spatial connectivity of 
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spins) and with Hamiltonians Ho and 'Hi as given in 
Eqs. (|4.2p and (|4.3p . The corresponding free energies 
are i^o,i = -^ In^^ cxp(-/3Ho,i); /3 = l/lfcsT) is the 
inverse thermal energy. We are interested in the free en- 
ergy difference AF = Fi — Fq. 

For the computation of this free energy difference we 
use the coupling parameter approach [55|. To this end 
we introduce the crossover Hamiltonian Her (A) which de- 
pends on the coupling parameter A, 



Hcr(A) = Ho + A (-Hi - Ho) = Ho + AAH, 



(4.4) 



with the Hamiltonian difference AH = Hi — Ho which 



in the 



,{h) 



Hf. 



The deriva- 



present context is Hn 
tive of the corresponding crossover free energy Fcr(A) = 
— 4 ln^^exp(— /3Hcr(A)) with respect to the coupling 
parameter reads 

, AHe-'3«"(^) 



dFer(A) _ 



dA 



^F^X)^ 



Ec 



Ec 



-/JW„(A) 



(AH)cr(A), 



(4.5) 

where (AH)cr(A) is the Hamiltonian difference AH aver- 
aged with respect to Her (A). Therefore one can compute 
the free energy difference by integrating (AH)^^ with re- 
spect to A: 



1 i 

AF = F,-Fo = J ^e'r(A)dA = J (AH), 



(A)dA. (4.6) 



For the forms of Ho and Hi given by Eqs. (|4!2|) and (j4?3| . 



respectively, one has 



AH = n\l"'^ - H^'^ , 



{h) 



,{h) 



-Her- = Wc7 = ^ST + (1 - A)HY^ + AHfi ^ . (4.7) 



2. Needle of broken bonds 

Here we consider the combination {i[0]j) correspond- 
ing to an "ordinary" needle which consists of an even 
number D of broken bonds in a strip with an even number 
W of rows. Figure [6] shows the example of a needle with 
D = 4 broken bonds in the center ofanLxW^ = 6x6 
strip. The Hamiltonians Ho and Hi for perpendicular 
and parallel needle orientation have the same form as 
the right hand side of Eq. (|4.ip but with reduced inter- 
action constants Ju.v;u'.v' which depend suitably on the 
coordinates {u,v) and {u' ,v') of nearest neighbor spins. 

Choosing without restriction um = 0, the needle with 
center at {Q,vn) is inserted in perpendicular orientation 
by "breaking", i.e., removing the D lattice bonds which 
at w = iiAr±l/2, ?;Ar±3/2, ..., vm^{D — 1)/2 extend from 
M = —1/2 to 7^ = 1/2. This is accomplished by adding 

D 

^± /^ = / , S-l/2,VN-(D+l)/2+k Sl/2,VN-{D+l)/2+k 

^ E (4-8) 

(inc.) 



A "I 



I I I I 



— Ju,v;u'.v> = 1 normal bond 
— bonds to boundaiies 
broken bonds not shown 
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(ir-i)/2 

3/2- 
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decreasing bond: 1 — A 
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(a) 



FIG. 6: Bond arrangements for a needle (O) of D = 4 bro- 
ken bonds with center (x) at {un,vn) = (0,0) at the midline 
u = of a strip with L = W — 6: (a) for the perpendic- 
ular orientation of the needle with Hamiltonian Ho; (b) for 
the parallel orientation of the needle with Hamiltonian Hi; 
(c) for the crossover Hamiltonian Her (A) which interpolates 
between (a) for A = and (b) for A = 1. Bonds of strength 
J — 1 are indicated by thin black lines. Broken bonds (J = 0) 
are not shown. The bonds indicated by green dashed and blue 
dot-dashed lines have strengths A and 1 — A which increase and 
decrease, respectively, as A increases. As in Fig. [T] the fluc- 
tuating spins Su,v = ±1 with \v\ < {W — l)/2 are indicated 
by empty circles while the fixed spins s„ ±(vt/+i)/2 = 1 in the 
two additional outside rows are indicated by full circles. The 
nearest neighbor bonds between fixed and fiuctuating spins 
are indicated by thick magenta lines. For the lower and up- 
per boundary i and j they have the strengths A^ and A ■ , 
respectively, where A'^' equals for an O boundary and ±1 
for a -\-/ — boundary. 



to the Hamiltonian Hst/i7 without needle (see 
Eq. (|4.2p ). Similarly for the parallel orientation of the 
needle one has to add 

D 

'^11 /^ ^ 2-^S_(D + l)/2+k,VN-l/2S-(D + l)/2+k,VN + l/2 
fc=l 

^ E (4.9) 

(deer.) 

so that AH = H^'^ - ufK Figures [^a) and (b) illus- 
trate these configurations for the special case vn = 0. 

The sums in Eqs. (|4.8p and (|4.9p have been charac- 
terized by subscripts (inc.) and (deer.) because in the 
crossover Hamiltonian following from Eq. (j4.7p . 



H(?)(A) = H(0) - A^ E -(1 - ^)^ E ' (4-10) 



(inc.) 



(deer.) 



they appear with a prefactor —A and —(1 — A), respec- 
tively, representing sums of products of spins coupled by 
nearest neighbor bonds with strengths XJ and (1 ~ X)J 
which increase and decrease, respectively, as A increases. 
Here H(0) =, Hst + J E +J E equals Hst in 

(inc.) (deer.) 

Eq. (|4.ip but with both types of nearest neighbor bonds 
missing which are broken in the perpendicular or the 
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parallel orientation of the needle. This corresponds to 
Fig. inijc) with both dashed and dash-dotted bonds re- 
moved (so that in 'w'-'^ only those bonds of Hst remain 
which arc outside a hole with the shape of a cross). In 

"Her (A), however, the two types of missing bonds are re- 
placed by the bonds of increasing and decreasing strength 
as illustrated in Fig. ^c) by the green dashed and blue 
dash-dotted lines. This obviously leads to the crossover 
from the perpendicular to the parallel needle orientation 
as A increases from to 1. 

On this basis, following the steps described by Eqs. 
(glll-dlTl) allows us to calculate F|[ -Fl = Fi -Fq = Ai^ 
for the combination {i[0]j). 



3. Needle of fixed spins 

In this subsection we consider the lattice version of the 
case (i[-F]j) in which a needle consisting of an odd num- 
ber D of spins fixed in the -I- direction is embedded in a 
strip with an odd number W of rows. Figures El^a) and 
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FIG. 7: Bond arrangements for a needle of D = 5 spins fixed 
in the -I- direction with the needle center (x) at {un,vn) = 
(0, 0) at the midline u = of a strip with L = 8 and W — 7 
for the perpendicular direction of the needle (a) and for the 
parallel direction of the needle (b). The strengths A'^' of the 
bonds near the strip boundaries are as explained in Fig. [S] 

(b) show the example of a needle with D — 5 in the center 
ofanLxVF = 8x7 strip. The partition functions and 
free energies of these lattice models remain unchanged 
if the bonds between the fixed needle spins and their, in 
all, 2D + 2 fluctuating nearest neighbors are removed and 
these 2D + 2 neighbors are coupled instead with the bulk 
strength ^J to a single exterior spin sg = +1 which is 
kept fixed in the + direction. The coupling to this single 
spin or to the D fixed spins has the same effect on the 
fiuctuating spins, namely that of a magnetic field acting 
on the 2D + 2 neighboring spins. Once this coupling to 
So is in place, for the following it is convenient to replace 
each of the D fixed needle spins by a freely fluctuating 
spin, i.e., free of any couplings. This changes the free 
energy per ksT only by _Dln2, independent of the orien- 
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FIG. 8: (a) and (b) describe how to mimic the configurations 
discussed in Fig. [7] by means of couplings (north-east arrows) 
to a fixed external spin so = -1-1; (c) shows the arrangement of 
bonds in the crossover Hamiltonian ?{„ (A) (see Eq. (|4.19|) ) 
with their strengths indicated in the right margin. The ar- 
rangement reduces to that of (a) and (b) for A = and A = 1, 
respectively. For further explanations see the main text. 



tation of the needle and thus drops out of AF. For the 
above example discussed in Fig. [7] this alternative model 
is illustrated in Figs.[5{a) and (b), with the couplings to 
the external spin denoted by north-east arrows. 

The corresponding additional terms in the Hamilto- 
nian can be written as 



U 



(+) 



and 



U 



(+) 



E+E-E-E 

(zero) (inc.) (one) (deer.) 



.._- {+) (+) 

E+E -E-E 

(zero) (deer.) (one) (inc.) 



(4.11) 



(4.12) 



/(+) 



(+)^ 



(so that AH = "H I ^ - "HY ) where the sums 

/ ^ = S0,iijv('S-l,iiiv+'Sl,-ujv+^0,iijv-(-l + '50,'Ujv-l)i (4.13) 
(zero) 
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(D-l)/2 



X! " X! [^0,VN+k{S-l,v,^+k + Si^y^+k) (4.14) 



(inc.) 



k=l 



and 



{D-l)/2 

^ = Yl [Sk,vt,{sk,VN+l+ Sk,vr,-l) (4.15) 

(deer.) k—1 

+ S-fc,t,jv(.S_A;/ujv + l + S-fe/ujv-l) 

~r^A:,t?Af ^^+l,^w ' ^ — k^VN ^~k — l,vj^ 



contain products of nearest neighbor lattice spins while 
the sums 

(+) 

(one) 



( + ) 

E 

(deer.) 



^0,VN + iD-l)/2-\-l + •^0,i;iv-(-D-l)/2-l 



"^"■5— l.lJiV "•" ^^,vn 



(4.17) 



(D-l)/2 
A:=2 



^ — l^vjsi+k ~T~ 5l,i;7v+fc ~r S_i^^jy_/^ + Si^yj^ — k 



and 



(+) 

/ , = ■5(D-l)/2+l.t)jv + *-(I5-l)/2-l.DjN 
(ine.) 



+ SO,i.N + l + ■SO.t.jv-l (4.18) 

Sfc/Ujv + l + Sfc,tlAr-l + S-k,VN + l + S-fc,tlAr-l 



(D-l)/2 

E ' 

fc=2 



contain products of a lattice spin and the fixed exter- 
nal spin So (which has the value sq = +1 and thus does 
not appear in Egs. (|4.16p - (|4.18p '). Here ^ contains the 

(zero) 

products of the center ( x ) spin of the needle with its four 
nearest neighbor spins which correspond to lattice bonds 
broken in both the perpendicular and parallel needle ori- 
entation (compare Figs. [SJa) and (b)). The products 

in ^ and ^ correspond to the remaining nearest 

(ine.) (deer.) 

neighbor bonds which are broken in the perpendicular 
and parallel needle orientation, respectively. The four 

(+) 
terms in ^ correspond to the bonds between the ex- 

(one) 

tcrnal spin and those four lattice spins which arc coupled 
to it in both the perpendicular and the parallel needle 
configuration (compare Figs. [Slja) and (b)). The sums 



(+) (+) 

^ and ^ contain the terms which correspond to 



(deer.) 



(ine.) 



the rest of the bonds to the external spin in the perpen- 
dicular and the parallel needle orientation, respectively. 
As in the previous subsection the notation for the vari- 
ous sums reflects the modulus of their corresponding pref- 
actors in the crossover Hamiltonian 

nipw - H(+) - Aj £ - (1 - A) J g 

(inc.) (deer.) 

(+) (+) (+) 

-^E-^^E-(i-^)'^ E • (4-19) 

(one) (inc.) (deer.) 

The first term n<-+^ = Ust + J E +JT,+JT^ 

(zero) (ine.) (deer.) 

corresponds to a strip with a cross-shaped hole where 
the bonds belonging to E j E i ^^'i E ^^^ miss- 

(zero) (ine.) (deer.) 

ing (i.e., in Fig. M,c) this means that the dashed and 
dash-dotted bonds and all the arrows are removed). In 
"Her (A) the contributions due to the two last types of 
bonds missing in 'H'"''', i.e., E ^^d E carry the pref- 

(ine.) (deer.) 

actors —A J' and —(1 — A)j7 of increasing and decreasing 
strengths , respectively. Moreover, bond contributions 
are added which couple the lattice spins contained in 

(+) (+) (+) 

E , El ^n-d E with strengths 1, A, and 1 — A, 

(one) (ine.) (deer.) 

respectively, to the external spin sq = 1. For L = 8, 
W = 7, D ~ 5, and the needle center x at wat = at 
the midline of the strip, the various bond strengths in 
Her (A) are shown in Fig. |8l^c) which clearly illustrates 
the crossover from the perpendicular to the parallel nee- 
dle orientation considered in Figs. [Slja) and MJ^) as A 
increases from to 1. 



4-. Details of the numerical implementation 

For the sequential generation of system configurations 
we have used the hybrid Monte Carlo method [53. One 
step consists of updating a Wolff cluster [5711 followed by 
L X 14^/4 attempts of Metropolis updates [53| of randomly 
chosen spins and of additional {D + 3)^ updates of ran- 
domly chosen spins in the square of size (D -I- 3) x (Z? -I- 3) 
with the center at position (0,vn). 

In order to determine the dependence of the free en- 
ergy on vn wc have used system sizes 1000 x 100 and 
1000 X 101 for (i[0]j) and {i[+]j) needles, respectively. 
For thermalization we have used 1.5 x 10^ MC steps, fol- 
lowed by the computation of the thermal average using 
8 X 10^ MC steps. These latter MC steps have been split 
into 16 intervals which facilitates to estimate the numer- 
ical inaccuracy. 

In order to determine the aspect ratio dependence of 
the free energies we have used various numbers of MC 
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steps (split into 8 intervals for estimating again the nu- 
merical inaccuracy) for various system sizes, varying from 
6 X 10^ MC steps for L = 4000 to 2.4 x lO^ for L = 200. 
We have used one fifth of these MC steps in order to 
achieve thermalization. 

Concerning the numerical integration over the 
crossover variable A, for every selected set of parame- 
ters (i.e., type of needle and boundary conditions {i[h]j), 
L, W, D, and vn) we have performed computations for 32 
points Afc = j^, fc = 0, 1, . . . , 31 and then we have car- 
ried out the numerical integration by using the extended 
version of Simpson's rule. 



V. COMPARISON OF ANALYTIC RESULTS 
WITH SIMULATION DATA 

Here we compare the "small needle" predictions from 
Sec. HI] for the quasi-torque AF with corresponding re- 
sults obtained by the Monte Carlo simulations described 
in Sec. IIVI On one hand this allows us to assess the 
performance of the truncated form of the "small needle 
approximation", i.e., to determine the smallness of the 
mesoscopic needle length D needed in order to be able 
to neglect higher order terms in this expansion. On the 
other hand, good agreement signals that all the meso- 
scopic distances and lengths, including D, chosen in the 
simulations turn out to be large enough for the lattice 
system to lie within the universal scaling region. 



A. Aspect ratio dependence for the double 
periodic strip 

For a needle in a double periodic strip, one can study 
the full dependence of AF (Eq. (|2.15p ) on the aspect 
ratio W/ L. Due to the symmetry of these boundary con- 
ditions, AF keeps its modulus but changes its sign as the 
values of W and L are exchanged. This implies that AF 
vanishes for W = L. Within the small needle approxima- 
tion, its quantitative behavior follows from the remarks 
below Eq. ([238]) yielding (see Eqs. (f2J6l) - (f238| ) 



AF 



AFi+AF^i D \ 

= "^ 2VK ^^(^/^)x 



keT 



1 n F> 



2 2 2W 



-JP{w/l) 



(5.1) 



for an {"ordinary", "normal"} needle with At ~ 
{1/2, -1/2} (see Eq. ([23])). Equation dSj]) is con- 
sistent with the sign change mentioned above because 
both $^^^ (see Eqs. ((Ml) and dH])) and (e)^^^ = 
W-^ff\w/L) (see Eq. dSH]) with x, = 1) re- 
main unchanged upon exchanging the values of W 
and L, implying Ap{W/L)/W'^ = ~Ap{L/W)/L'^ and 
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FIG. 9: Normalized quasi-torque AF = F|| -Fj_ (Eq. ^^) 
acting on a small mesoscopic needle of length [23] 13 in a dou- 
ble periodic strip P as a function of the aspect ratio W/L 
of the strip. AF > implies that there is a preference for 
the perpendicular orientation. AF vanishes at W/L = 1 and 
is antisymmetric around this point. The black dotted line 
denotes its limit for W/L = 0. For "ordinary" [O] and "nor- 
mal" [-I-] needles with 13 = 20 and D = 21, respectively, the 
simulation data (green squares and magenta diamonds) agree 
quite well with the analytic prediction (red dash-dotted line) 
according to the first term on the rhs of Eq. (|5.1|l . The fact 
that the (green) squares lie above the (magenta) diamonds 
agrees with the tendency arising from the second term AF„; 
f{P) 



in Eq. (|5.1fl due to /l < 0. However, the splitting between 
[O] and [-I-] due to AF„; is quite small. (Considering 13 = 20 
and W — 100, the explicit form of /i ' given in Eqs. (|A2|I - 
(|X4)) leads at, e.g., W/L = 0.2 and W/L = 0.5 to the values 
1 ± 0.003 and 1 ± 0.027, respectively, of the square bracket in 
Eq. ([STT})- A quantitatively reliable check of this small split- 
ting effect requires accurate data from lattice systems with 
larger values of 73, W , and L which are closer to the asymp- 
totic scaling limit than the ones presently available. For more 
details see the main text. 



fP{W/L)/W = /i 
pendix lAll — sgn(F - 



\l/W)/L. Since, due to Ap- 
W) X Ap{W/L) is positive and 



/i is negative, Eq. (|5.1|) implies the following: 

(i) Both "ordinary" and "normal" needles prefer to 
align perpendicular to the longer axis of the dou- 
ble periodic strip. (Concerning related effects see 
Ref. I!].) 

(ii) While in leading order oc D^ the strength of this 
preference is independent of the needle type /i, it is 
strengthened (weakened) for an "ordinary" ("nor- 
mal" ) needle by the correction of order D^ . 

A quantitative comparison of the aspect ratio depen- 
dence predicted in Eq. (|5.1|) with our simulation data is 
provided in Fig. ID This figure shows data [601 for the nor- 
malized quantity [AF / {kBT)]/ {D /WY within the range 
< D/W < 1/2. For the "ordinary" needle of "bro- 
ken bonds" the data were obtained from systems with 
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W = 100 and D = 10 and 20. The data for the "normal" 
needle of "fixed spins" stem from systems with W = 101 
and D = 11 as wcU as £) = 21. The data for D = 20 
and D = 21 are in rather good agreement with the ex- 
pression [AFi/{kBT)]/{D/W)^, as predicted according 
to Eq. (|5.ip and Appendix lAll for the leading order con- 
tribution ex D^, shown as the red dash-dotted hne. The 
data for the "ordinary" needle are indeed shghtly larger 
than those for the "normal" needle, as predicted by the 
next-to- leading order contribution ex D^ in Eq. (|5.ip . 
The stronger deviations of the data for D = 10 and 
D = 11 from the analytic approximation presumably in- 
dicate that these smaller needle lengths lie outside the 
mesoscopic scaling region required j24| for the validity of 
Eq. dSH). 
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B. Dependence of the free energy anisotropy on 
the spatial position of needles in strips. 

In strips with boundaries AF depends on the position 
vn of the needle in the strip. Here we consider strips of 
infinite length L, either with arbitrary boundaries {i-,j) 
and containing an "ordinary" needle or with boundaries 
(O, O) and containing a "normal" needle. In all these 
cases AF ~ AF; + AF„; is predicted to have the form 
given by Eqs. ([236)) - ([235)1 . 

Ordinary needles 

While AFi does not depend on vn and is given by 
Eq. (|2.17l) with A^.j given below Eq. (|2.9p . in the cases 
(i[0]j) one has 



AFnl 

knT 




(5.2) 



which depends on w^v via the following simple expressions 
for ffij! 

go^o = 3(cosy)"3 - (5/3)(cosF)-i 
.9+,+ = -9o,o 

9+.0 = -9o,+ = (tanF)[3(cosV^)-2 + (1/3)] 
5+,-=5-,+ = -[3(cosy)-3 + 

+ (7/3)((cosF)-^-4cosy)] (5.3) 

with V = ttvn /W. They follow upon inserting /e from 



A^ 



(O) 



1/2 (see Eq. ([23|)) into 



Eq. ([AT4I1 and Ai'''' 
Eq. ^M- 

Figure [lU] compares our simulation data for the cases 
(*[^]j) with the corresponding analytic predictions for 
various positions vn of the needle and for various bound- 
ary conditions {i,j) of the strip. The plots show the data 
for [{AF-AFi)/{kBT)]/{D/W)^ with AF obtained from 
a system with L = 1000, W = 100, and D = 20. The 
comparison with [AFni/{kBT)]/{D/Wy = c,;.j following 
from Eqs. ()5.2I) and ()5.3p and shown by lines is very fa- 
vorable (i.e., AF is captured well by AF; -I- AF„i), except 



FIG. 10: Quasi-torque AF acting on a small mesoscopic "or- 
dinary" needle O in a long strip with boundaries (J, j) cor- 
responding to the cases {i[0]j). This plot shows its depen- 
dence on the position vn of the needle in the strip and on the 
boundary conditions {i,j)- The lines are suitably normalized 
expressions Cij (see Eqs. (|5.2p and (|5.3|l and the text below 
the latter one) for AF„; which, according to Eq. (|5.3|) . are 
symmetric around vn/W — for {i,j) = {0,0), {+,+), and 
(— ,-f) but antisymmetric for {i,j) — {0,+). In order to ob- 
tain the numerical data we have computed AF by means of 
the coupling parameter approach (Sec. lIV|) and then have sub- 
tracted the analytic expression (|2.17|l for AF;. There is very 
favorable agreement of the simulation data, which correspond 
to W/L = 0.1, with the analytic predictions for W/L = 0. 
The case (— [0]-f) is exceptional in that for it the agreement 
is only fair, i.e., in this case there are sizeable corrections to 
AF — AFi + AF„i beyond AF„; . The statistical error bars 
are comparable with the symbol sizes and tlierefore they are 
omitted. 



for the case (—[0]+) in which it is only fair. According to 



-ij 



(7r3/256)(D/IF)g. 



«j- 



Eqs. ([Ol and ([O)) one has 

In order to visualize the approach of the limit of in- 
finite strip length F = 00, Figs. [11] and [T^] show [60| 
the dependence of the quasi-torque on the aspect ra- 
tio W/L for needles of length £» = 10 or F = 20 with 
their center v^^ = at the midpoint of a strip of width 
W := 100. The cases {0[0]+) and {0[0]~) should lead 
to the same AF, due to the (+ ^ — ) symmetry [3J| of 
the Ising model, and AF„; should vanish for vn = 0, due 
to Eq. ()2.18p and Ref . [H] . These properties are reflected 
rather well by the data in Fig. [TT] for D = 20 (shown as 
squares and diamonds) which are close to the values of 
[AFi/{kBT)]/{D/W)^, equal to -7rV96 for W/L = 
(shown as the horizontal straight line) and about 5% 
smaller for W/L ~ 1/2. Here we discard the data for 
W/L < 0.1 which carry large error bars. This should 
be compared with the case {0[0]0) discussed in Fig. [T^ 
where in a strip of width W ^ 100 the data for the nee- 
dle of length D = 20 (shown by squares) attain much 
closer the asymptotic value (7r2/192)(l -|- 207r/100) of 
[(AF; + AFni)/{kBT)]/iD/W)^ predicted for vn ^ 
and L = 00 (shown as the uppermost straight line) as L 
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FIG. 11: Dependence of the quasi-torque AF on the aspect 
ratio W/L of the strip for iO[0]+) and iO[0]-) bound- 
ary conditions with the needle at tlie strip center vm ~ 0. 
In these cases the next-to-leading contribution AFni van- 
ishes [431 for arbitrary W/L and the simulation data for 
D = 20 (squares and diamonds) are indeed close to the val- 
ues given by —ttAo,+{W/L)/4: (full black line) correspond- 
ing to the leading contribution AFi which are — 7r^/96 for 
W/L = and about 5 percent smaller for W/L — 1/2 (see 
Eqs. (|A5|) . (|A7|) . and (|A9|) '). This decrease is too weak to be 
reflected by the present simulation data. Due to the inher- 
ent {+ O -) symmetry, one has AFf°[°l+' = AF''°^°^-'> for 
arbitrary values of D, vn, W , and L. This exact identity is 
embodied in the form of the corresponding lattice Hamiltoni- 
ans AT-L and "Her (A) in Sec. IIVI However, the results for the 
thermal average {A?^)ci(A) and its integral AF in Eq. (|4.6|) . 
calculated by means of the statistical Monte Carlo method, 
violate this identity within the numerical inaccuracy. Note 
that the ensuing deviations in the above data are of tolerable 
size, at least for the larger value of D. 



increases, i.e., as W/L becomes smaller. 

Normal needles in strips without broken symmetry 
Now we consider the (equivalent) cases (0[-|-]0) and 
(0[— ]0) of a "normal" needle with ft. = -I- or — 
in an (O, O) strip. Figure [13] shows data [601 f^^' 
[AF/(kBT)]/{D/Wf for such a needle of length D = 21 
at various positions wjv in a strip with L = 1000 and W ~ 
101. They compare favorably with the corresponding an- 
alytic expression -(7r/4)Ao,o + (7i"^/256)(i:'/Ty)(-go,o) 
for [{AFi + AFni)/{kBT)]/{D/Wf with L = oo, which 
is shown as dashed line. This expression follows from 
Eqs. (|2.16p - (|2.18p by repeating analogously the line of 
arguments leading to Eqs. (|5.2I) and (|5.3p . Here the sign 
in front oi go,o differs from that in Eq. (|5.2p for the "ordi- 



nary" needle because A 
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AP (Eq. ([H 



In strips of infinite length L the pseudo-torques AF{vn) 
for the present case (0[-|-]0) and for the case (+[0]+) 
considered above are identical. This follows from a dual- 
ity argument similar to the one in Ref. [4J], is in agree- 
ment with the expression AF; + AFni in Eqs. (j2.16p - 
(|2.18|) . and is quite well reflected by the simulation data 
in Figs. [TU] and [T31 Concerning the aspect ratio depen- 
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FIG. 12: Similar as Fig. [TT] but for the cases {0[0]0) and 
{0[+]0) in which, forL = oo, [{AFi+AF^i)/{kBT)]/{D/Wf 
equals {n^ /192)[1 + ttD/W] and {n^ /192)[1-ttD/W], respec- 
tively, so that the leading contributions (dashed line) are the 
same. For small W/L the simulation data for D = 20, W — 
100 and D = 21, W = 101 (green squares and magenta dia- 
monds) are indeed close to the corresponding uppermost and 
lowest horizontal line, respectively. 
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FIG. 13: Dependence of the quasi-torque AF on the position 
vm of the needle in the (0[+]0) case. The simulation data 
(triangles) agree well with the analytic prediction (dashed 
line) given by AFi + AF^i (see the main text). 



deuce, for W/L < 0.1 the data for the case {0[+]0) of 
a "normal" needle with vn ~ and D = 21 (shown by 
diamonds in Fig. 112^ have the tendency to approach the 
predicted hmiting value (7r^/192)(l — 2l7r/101) (indicated 
in Fig.[T2|by the lowest horizontal straight line). 



C. "Normal" needles in strips wth broken 
symmetry 

In this subsection we consider a "normal" needle with 
h = + embedded in strips with at least one "normal" 
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FIG. 14: Quasi-torque AF acting on a + needle in tlie cen- 
ter of a ( — ,— ) strip. We show simulation data for needles 
of length D = 11 and 21 in strips of width W — 101 and of 
various lengths L. The data for D = 21 (diamonds) agree 
quite well with the "sna" prediction —0.4143 (dash-dotted 
line) for a strip of infinite length as explained in the first two 
paragraphs of Subsec. IV CI For the smaller needle of length 
D = 11 there is a significant discrepancy between the simula- 
tion data and the "sna" prediction -0.3072 (dashed line) for 
a strip of infinite length. 



boundary. These are the boundary conditions (i,j) = 
[{+, +), (-, -), (-1-, -), (+, O), (-, O)] for which the order 
parameter profiles ((/))st = W-'^/^fj^^'^^ in Eq. dgT]) are 
nonvanishing and Q in Eq. (|2.13p contains a contribution 
ex D^'^. In these cases we do not compare our simulation 
data with the expanded analytic form of AF, because 
confining this expansion to the two leading powers of D, 
as in Eq. (|2.16l) , is expected to be insufficient for reaching 
agreement with presently accessible simulation data. As 
explained there, we rather set out to compare the simu- 
lation data with the full expression "sna" of the "small 
needle approximation" following from Eqs. (|2.10p - (|2.15p . 
While we disregard the terms of order D'^ denoted by 
the ellipses [4l[ in Eqs. (PTT^ and dHH]), we leave the 
logarithm in Eq. (j2.15p unexpanded. 

We illustrate this point for the case (—[+]—), i.e., a -I- 
needle embedded in a (— ,— ) strip. Figure [T4l shows [6G| 
that [AF/{kBT)]/{D/W)'^ as obtained from simulations 
for a needle (with D = 21 and located in the center of 
a strip with W = 101) agrees, for smaller aspect ratios 
W/L, quite well with the prediction —0.4143 for L = 00 
(see the diamonds and the dash-dotted bottom line in 
Fig. I14p . The prediction follows by inserting the present 
value e ^ 2l7r/101 == 0.653 of 6I ee ttD/W into the "sna" 
expression AF/{kBT) = — ln(Z|j/Zj_) where 

Z||,^= 1 -(20)i/«-2-'^+i/«0i^/« + 2-30 + 2-^03 
+ (1,-1) X [-2-^(1/3)02 

-f 2-6+1/8(1/3)017/8 _ 2-i"303] , (5.4) 
see Eqs. (^3)1 and (pi^ - (PTi)) with the value A_ _ = 



-7r/48 taken from below Eq. ((2J|) as wefi as Eqs. (|A14p 
and (|A15[) . For this value of 9 the validity of the "sna" 
is confirmed by Fig. 21 as far as the contribution F|| to 
AF is concerned (compare the discussion in the para- 
graph below Eq. (|3.15p ). This should be contrasted with 
the prediction —0.0017 (as compared with —0.4143, see 
above) which would follow from AF being truncated af- 
ter the leading and the next-to-leading order in D, i.e., 
from (AF, + AF^O/keT = [{26)^ - {29)^^^]/ {3 x 2^) (by 
using again Eq. (j5.4p ) and which would disagree with the 
simulation data by a factor of about 250. 

The reasonable agreement between the simulation data 
for D = 21 and the "sna" persists for the dependence on 
the needle position vn shown as the squares and the green 
dashed line in Fig. 1151 provided |fAr|/VF remains small. 

However, for |w7v|/W^ ^ 0.15 an even qualitative de- 
viation develops in that the analytic approximation pre- 
dicts a point of inflection and a minimum which is not 
supported by the simulation data. We attribute this fail- 
ure to inadequacies of this approximation near the strip 
boundary which have been discussed in detail for F± in 
the half plane (see the paragraph in Sec. IIIII addressing 
Fig. El^a)). In particular the unphysical minimum of the 
green dashed line in Fig. [15] is related to the maximum 
of the "sna" result shown in Fig. El^a). 

This dependence on vn/W is shown in Fig. [TS] also for 
the other cases (-[+]0), (+[+]0), (+[+]-), and (+[+]+) 
of a + needle embedded in strips with at least one "nor- 
mal" boundary. Also in these cases the dependences of 
the data on vn/W are reproduced reasonably well by the 



sna 



VI. SUMMARY AND CONCLUDING 
REMARKS 

A critical solvent such as a binary liquid mixture at 
its continuous demixing transition induces a long-ranged, 
so-called critical Casimir interaction between immersed 
particles which is universal on mesoscopic length scales. 
For nonspherical particles the interaction depends not 
only on their sizes and distances among each other but 
also on their shapes and mutual orientations. We have 
studied a critical system belonging to the Ising bulk uni- 
versality class in two spatial dimensions and particles of 
needle shape. 

As described in the Introduction we have considered 
Ising strips at the bulk critical point with an embedded 
needle which is oriented either parallel or perpendicular 
to the symmetry axis of the strips. For such systems 
our analysis benefits from the wealth of knowledge ac- 
cumulated for the two-dimensional Ising model at criti- 
cality and the comparative ease to implement a needle 
along a row or column of the square lattice. The cor- 
responding effective interaction between the needle and 
the two strip boundaries is probably the simplest exam- 
ple to check by Monte Carlo simulations the range of va- 
lidity of asymptotic analytic predictions for the orienta- 
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tion dependence of critical Casimir interactions between 
nonsplierical mesoscopic particles which go beyond the 
Derjaguin approximation. 

In Sec. HV] and Figs. [H El [3 and [8] we explain how to 
implement in the lattice model boundary properties at 
the two confining surfaces of the strip and for the needle, 
which locally induce one of the two demixing bulk phases 
("normal" surface universality classes + (or — )) or in- 
duce disorder and suppress demixing ( "ordinary" surface 
class O). The geometrical features of the corresponding 
continuum description are explained in Fig. [2] 

Primarily, we analyze the free energy AF = Fu — F±_ 
required at bulk criticality to turn the needle from an 
alignment perpendicular to the strip (with free energy 
F±) to parallel alignment (with free energy F||). Thus 
one has AF > and AF < if the needle prefers the 
perpendicular and parallel alignment, respectively. AF 
depends on the length D of the needle, the width W and 
length L of the strip, the distance vn of the needle center 
from the midline of the strip, and the surface universality 
classes h of the needle and {i,j) of the two boundaries of 
the strip. Here we consider needles of small mesoscopic 
length D for which predictions are available from the so- 
called "small needle expansion" explained in Sec. [Ill 

Before presenting below an itemized summary of the 
quantitative comparison with the simulations we point 
out a few qualitative observations. For the needle it is 
advantageous to reside in a spatial region and take an 
orientation which suits its boundary condition. 

(i) First, consider a needle of universality class h in a 
half plane with the boundary belonging to surface univer- 
sality class i. The needle will prefer the vicinity of this 
boundary if i = ft,. Since in the case i = + (i = O) 
the boundary-induced order (disorder) increases upon 
approaching the boundary - as described, according to 
Eq. (|2.5p . by the concomitant increasing density profiles 
of the order parameter (of the energy) in the half plane 
without needle - this leads to an attractive force between 
needle and boundary. Since both increases are stronger 
than linear, for a fixed needle center in both cases the nee- 
dle will adopt an orientation perpendicular to the bound- 
ary of the same class. Likewise, for different universality 
classes i ^ h oi the needle and the boundary the force 
will be repulsive and for a fixed center the needle will 
orient parallel to the boundary [lj|. These orientational 
preferences depend on the needle position, increasing and 
decreasing, respectively, with decreasing and increasing 
distance ajv of the needle center from the boundary. 

(ii) In a strip there is an additional qualitative effect 
in that there is a contribution to the orientational pref- 
erence of the needle which is independent of its position 
Vn in the strip but depends on the combination {i,j) 
of surface universality classes of the two strip bound- 
aries. In a long strip this contribution favors an align- 
ment of the needle perpendicular (parallel) to the bound- 
aries if they belong to the same (different) universality 
class i = j (i 7^ j). For the double periodic strip it 
favors an alignment perpendicular to the longer axis. 



This is reminiscent of - and related to - a well known 
corresponding anisotropy [59| of the two-point averages 
{0{r - s/2)0(r + s/2))st in the strip of the densities 
O = ^ of the order parameter and O = e of the energy. 
For small mesoscopic distances |s| the anisotropy causes 
the two-point averages to be larger for s perpendicular 
(parallel) to the two strip boundaries ii i = j (i ^ j) 
and for s perpendicular to the longer axis of the double 
periodic strip. 

All of these expectations arc confirmed and further re- 
fined by the quantitative predictions of the "small needle 
expansion" . It predicts, in particular, that the anisotropy 
AF due to [61| the aforementioned effect (ii) depends, 
apart from the boundary classes (i, j) of the strip, only 
on the length and not [40| on the universality class 
(/i = +, — , or O) of the small mesoscopic needle. 

In the following three blocks A, B, and C we list our 
main results. 

(A) We have demonstrated that the quantitative pre- 
dictions of Sec. inifor the rich structure of the orientation- 
dependent interaction of non-spherical particles with a 
"small mesoscopic size" can actually be observed in a lat- 
tice model. This is a nontrivial result in view of a twofold 
size condition: the particle size being small compared 
with other geometric features and being large on the scale 
of the lattice constant. Our Monte Carlo simulations 
with needle lengths of about 20 lattice constants [23] have 
the potential to closely approach the asymptotic regime 
of the "small mesoscopic needle" , leading to results for 
AF which agree quite well with the analytic predictions, 
without adjusting any parameter. At the same time, the 
parameter range in which there is good agreement pro- 
vides a (conservative) estimate for the range of validity 
of the "small needle approximation" ("sna"), which is a 
truncated form of the systematic expansion for the uni- 
versal quasi-torque AF in terms of the needle size. 

Al. For a strip with double periodic boundary condi- 
tions, for which the aforementioned effect (i) is ab- 
sent, our data for AF reproduce, for the larger nee- 
dle lengths D — 20, 21, quite well the predicted 
dependence of AF on the aspect ratio W/L of the 
strip and the predicted independence of AF of the 
universality classes h = O and h — + oi the needle, 
in leading order (see Fig. [9] and the discussion in 
SubsecEH). 

A2. In strips with actual boundaries, effects of both 
type (i) and (ii) are present and the former ones 
lead to a dependence of Fy, F±, and AF on the 
position Vn of the needle center within the strip. 
Here the analytic predictions hold for strips of in- 
finite length L ~ oo and Figs. [Til [T^l and [HI show 
the dependence of the simulation data on the as- 
pect ratio W/L of the strip. 

A2a. For the five possible and relevant combinations 

i^[h]J) = iO[0]0), i+[0]+), (-[0]+), iO[0]+), 
and (0[-|~]0) of universality classes of the two 
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FIG. 15: Quasi-torque A_F acting on a small niesoscopic "nor- 
mal" needle + in a long strip in which one or both bound- 
aries break the Ising symmetry. We show the dependences 
of A_F on the position vn of the needle in the strip and on 
the boundary conditions (i, j) of the strip. There is reason- 
able agreement of the simulation data for W jL = 0.1 with 
the analytic predictions of the "small needle approximation" 
for W jL — which are discussed in the first paragraph of 
Subsec. IVCI For the case ( — [+] — ) the "sna" can be trusted 
only for |wjv |/W^ smaller than « 0.15 (see the discussion in the 
last but one paragraph of Subsec. |VC|) . The statistical error 
of the simulation data is comparable with the symbol sizes. 



boundaries of the strip and of the embedded nee- 
dle, for which the expansion in Eqs. (|2.16p - (|2.18p 
of AF in terms of powers of D is appropriate, the 
dependence on wjv of our data for AF is discussed 
in Subsec. rVBI In Figs. [TOl and fT3l this dependence 
is compared with the approximate analytic predic- 
tions. There is fair agreement for the case (— [0]+) 
while in the four other cases the agreement is very 
good. 

In the case (0[0]+), due to effect (i) one has 
AF{vn) > and AF{vn) < near the O and + 
boundary, respectively, so that the needle prefers 
the perpendicular and parallel orientation, respec- 
tively. Since effect (ii) for unequal strip boundaries 
turns the needle parallel to the strip, the range of 
vn values allowing for parallel orientation is the 
larger part of the accessible values of vjy, i.e., the 
point Vn where AF changes sign is closer to the O 
boundary. 

In the strip of infinite length L ~ oo duality ar- 
guments [3^|4J] predict for the cases (+[0]+) and 
(0[-|-]0) the same expression for AF{vn) which 
within the "sna" (Eqs. (PTB)) - (PT5)) ) is visualized 
as the full line in Fig. [TU] and as the dashed line 
in Fig. [131 It predicts that in this case the nee- 
dle prefers the orientation perpendicular to the 
strip for small |wjv| where effect (ii) prevails so 
that AF > 0, while for larger \vi\[\ where effect 



(i) prevails it prefers the parallel orientation with 
AF < 0. The good agreement with the data in 
Figs, [in] and [T3| tells again that within our lattice 
model one can access the scaling region and thus 
capture the corresponding universal small particle 
behavior. 

A2b. Unlike the combinations considered above 
(in A2a.) the five combinations {i[h]j) ~ 

(+[+]+), (-[+]-), (+[+]-), (+[+]o), {-[+]0) 

discussed in Subsec. IV CI involve a needle which at 
Tc is subject to the order parameter profile induced 
by the strip boundaries {i,j) via its nonvanishing 
amplitude A^^^ (see Eqs. ([^m)) - (|^Ti)) ). This leads 
to a contribution to the partition function ex D^'^ . 
For the comparison with the simulation data for 
AF, in these five cases the corresponding "small 
needle approximation" should not be implemented 
by directly expanding AF in terms of powers of D 
but by expanding the corresponding partition func- 
tions (see the remarks above Eq. (|2.16p and in the 
first paragraph of Subsec. [VCj. The comparison of 
the analytical predictions for these cases with the 
corresponding MC data is shown in Fig. [T5] The 
orientation preferred by the needle (i.e., the sign 
of AF) is in conformance with [62| the qualitative 
observation (i); also quantitatively the MC data 
agree quite well with the predictions. According 
to Fig. [151 in the case (—[+]—), however, the good 
agreement for small values of |wAr|/VF does not 
extend beyond It^Arl/VF « 0.15. In this case the 
expected breakdown of the "small needle approxi- 
mation" near strip boundaries starts already rather 
close to the strip center producing an unphysical 
minimum of AF{vn) at |wAr|/W^ ~ 0.25. 

(B) We have investigated the effective interaction be- 
tween the needle and the boundaries also for needles of 
arbitrary mesoscopic length D, i.e., beyond the regime of 
the "small needle expansion" . 

Bl. These exact results have been derived in Sec. IIIII 
and Appendix [Bl for the free energy of a needle in 
the half plane with an orientation perpendicular to 
its boundary i and of a needle embedded along the 
midline of a strip of infinite length and with the 
two boundaries being members of the same sur- 
face universality class [i ~ j). These geometries 
are related to the W x L strip without a needle 
(Appendix IA2p via conformal transformations of 
the Schwarz-Christoffel type [63|. For various com- 
binations (/i, i) of the surface universality classes 
these effective interactions are shown in Figs. [3J51 
and display the crossover between the small needle 
regime and the limiting behavior for which the nee- 
dle approaches and nearly touches the boundary or 
becomes much longer than the width of the strip. 

B2. Apart from the importance in their own right these 
universal results allow us to better understand the 
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limitations of the "small needle approximation". 
For example, in Fig. [Sfa), for {i[h]) = (—[+]) 
this approximation reproduces the exact result for 



F 



(-[+]) 



very well down to a distance of the closer 



end of the perpendicular needle from the boundary 
corresponding to half its length (a< = D/2 = aiv/2 
so that 1? = D/{2a]y) = 0.5). For even smaller dis- 
tances the approximation develops an unphysical 
maximum, which is the counterpart of the mini- 
mum produced by the approximation for AF in the 
case (—[+]—) (see Fig. [13]). In contrast, the exact 
result crosses over to a logarithmic increase which 
diverges when the needle "touches" the bound- 
ary. We recall that the validity of the "small nee- 
dle approximation" requires that the mesoscopic 
length D of the needle is small compared with 
the distance a< between the closer end of the nee- 
dle and the boundary, whereas the exact result is 
valid if only the microscopic lengths are sufficiently 
small compared with a<. Concerning the free en- 
ergy i^ii ' '"' for embedding a + needle extending 
along the midline of a (— , — ) strip, Fig.|l]exphcitly 
demonstrates the ensuing improvement if one ex- 
pands the corresponding partition function rather 
than the free energy in terms of the needle length 
(see the remarks above Eq. (|2.16l) and in the first 
and second paragraph of Subsec. IV C[) . 

(C) Our results allow us to conclude which of the fea- 
tures studied here are of a more general character and 
thus can be expected to show up also in spatial dimen- 
sion d = 3 and which ones are specific for d — 2. 

CI. In dimensions d > 2 nonsphcrical mesoscopic par- 
ticles such as uniaxial ellipsoids [l^, [iJ] or dumb- 
bells of two touching or interpenetrating spheres of 
equal size [13| in a half space and interacting with 
its planar boundary wall have been considered for 
various universality classes {i[h]). Both for particle 
sizes small [ij, [6J] and comparable [Ijl with re- 
spect to the distance from the wall, these particles 
prefer, at fixed particle center, orientations perpen- 
dicular (parallel) to the boundary of the half space 
a i ^ h {i y^ h). These are the same preferences as 
described in paragraph (i) above for our small nee- 
dles in the half plane. In the mean-field treatment 
of Ref. [ij] it was pointed out that at fixed closest 
surface-to-surface distance of particle and wall the 
preferred orientations display the opposite trend, 
i.e., being parallel (perpendicular) to the bound- 
ary of the half space ii i = h {i y^ h). These lat- 
ter preferences are in agreement with our finding 



i^lh]) 



F 



(t[h]) 



in d ~ 2 that the free energy F, 
quired to turn a long needle in the half plane about 
its closer end from the perpendicular to the par- 
allel orientation, is dominated by the first term, 
fiP''*'^ = kBTA,j,D/a< (note the signs of A,,,, 



second term, given by Eq. p.lOp . depends only log- 
arithmically on the large ratio £'/a< of the needle 
length D and the distance a< of the closer end of 
the needle from the boundary and thus can be ne- 
glected relative to the linear increase exhibited by 



m]) 



II 

given below Eq. 



The reason is that the 
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C2. However, there are also effects which are specific to 
two dimensions, due to the symmetries based on the 
duality transformation and the much wider class of 
conformal mappings. For example, the aforemen- 
tioned equality of the particle insertion free ener- 
gies as function of vn for the cases (0[-|-]0) and 
(+[0]+) with infinite strip length - based on the 
duality symmetry [3^, |4J] of the d = 2 Ising model 
at the bulk critical point - has no correspondence 
for a particle between parallel walls in d > 2. Like- 
wise, the independence of h of the leading small 
particle contribution AFi to the quasi-torque AF 
in the case {0\h]0) (see Eq. ((2Tfl) ) is valid in d = 2 
but not in d > 2. The derivation of this result in 
d — 2 hinges on using a conformal transformation 
as described in Ref. [40] which is not available for 
a non-spherical particle in d > 2. 

C3. There is an interesting difference between "ordi- 
nary" needles in d = 2 and d = 3 which arises 
from the dependence on their width W. Here we 
have considered needles with a mesoscopic width 
W which is much smaller than the needle length 
D. For both "ordinary" and "normal" needles in 
d = 2 the effects they induce in the embedding 
critical system (such as the density profiles), and 
thus their effective interactions with the bound- 
aries leading to force and torque, only depend on 
their length D but not ^ on their width W. This 
applies also [6a] to a "normal" needle in d = 3. 
However, for an "ordinary" needle in d = 3 with 
fixed length D the strength of these effects de- 
creases [631 upon decreasing the width W. For ex- 
ample, within the small needle expansion the pref- 
actor of the energy density O = e is proportional 
to DW'-^ = DW°-^^ in d = 3 while in d = 2 it 
is proportional to D-^' — D and independent of W 
(see Eq. ([IJ) and Ref. (H). 

Extending our detailed investigations to non-spherical 
particles in three spatial dimensions and for the whole 
neighborhood [6g of the critical point is desirable but 
beyond the scope of the present study. For Monte Carlo 
simulations in d = 3, studies of (square shaped) disks 
and of "normal" needles with their properties being in- 
dependent of their width look most promising |65| . Cor- 
responding quantitative analytical predictions remain a 
challenge. The "small particle expansion" can in prin- 
ciple be extended to, e.g., a circular disk. However, the 
amplitudes and averages of the corresponding operators 
can be obtained only approximately. A mean field treat- 
ment corresponding to d = 4 is given in Ref. |13j but al- 
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ready one-loop field theoretic calculations, corresponding 
to the first order contribution in an expansion in terms 
of 4 — d, look quite demanding. 
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Appendix A: STRIP WITHOUT NEEDLE 

1. Double periodic boundary conditions 

The behavior near bulk criticality of finite Ising strips 
with periodic boundary conditions in both Cartesian di- 
rections has found a long lasting theoretical interest, 
starting to the best of our knowledge in 1969 with the 
seminal paper [46[ by Ferdinand and Fisher. Therein 
the aspect ratio dependence of the free energy is given 
explicitly and was re-derived by other methods in later 
studies. The quantity — $sx addressed below Eqs. (|2.18p 
and (|5.ip is given by the scale free part of In Z^m. h^ 
Eq. (3.37) of Ref. Hg. It is equal to InZ^ in Ref. [13 
and to In Zpp for the Ising model in Ref. [43]. Adopt- 
ing Cardy's notation (sec Eqs. (2.4)-(2.6) and Table 1 in 
Ref. [431 ) the relation for the Casimir (or stress tensor) 
amplitude Ap given below our Eq. (|2.18p yields the result 

ApiW/L) = Ap{l/S) = -7T/12^{d/dS)ln{[xiiiS)]' 

+ [X2liS)? + [X22i6)?} (Al) 

where for the Ising model one has the conformal charge 
c = 1/2 and the parameter m = 3 in the functions Xpq 
introduced by Cardy. 

The energy density (e)sT = fe {W/L)/W defined be- 
low Eqs. (|2.3[) and (|5.ip is of course independent of u and 
V and turns out to vanish |46l . |70[ for strips of infinite 

length L, i.e., /i (0) = 0. The asp ect ratio dependence 
has been determined in Refs. [4g and [43 and in our 
notation reads 

fP (W/L) EE fP {l/5) = -nE{5)/Z{S) (A2) 

where 

E{S) = [C/i/24n^^^(l - C/")]2 (A3) 

with U = exp(— 27r(5) and 

Z{5) =. {1/2){2E{25)/E{5) + [E{6)f /[E{25)E{6/2)] 

+E{5I2)IE{5)}. (A4) 



fiP) 



Our quantities /e , 5, E, and Z correspond to the 
quantities — (e)^, — JT, jry^j, and Z^ , respectively, in 
Ref. [431 (see the introductory remarks in Sec. IIIII and 
Eqs. ([S3, ([SH), and ([XTi)) therein). 



2. Strip with boundaries (i, j) 

For our strip (ST) with the aspect ratio 1/(5 = W/L 
and boundaries 



(z,j) = [(0,0),(+,+),(+,-),(+,0)] 



(A5) 



Cardy's results |43| for the partition functions and our 
Eqs. ([13 and ([12D yield 

A,,,(l/<5) = -V48 - (d/dS) Hxn{S/2) + X2i{S/2), 
XiiiS/2), X2i{S/2), X22{S/2)] (A6) 

for the universal stress tensor am plit udes in terms of the 
functions Xpq introduced in Ref. [ij] . For later reference 
we rewrite Cardy's expressions in the form 

A,j (1/(5) = (7r/48) [ -1 + 12K{d/dK) ln{ [Sn + S21, 

En, S21, E22]n,^^i(l - «4")-i} j (A7) 
and alternatively as 

A, J (1/(5) = (7r(5- Vl2) I 1 - 24a{d/da) ln| (l + 



+ [(5ii + 52i)/2, 5n, 521, 522]) fj (1 - ^r^")-!} j 

(A8) 



n=l 



which converge rapidly for long strips with large S (cor- 
responding to extended, closely spaced boundaries) and 
for short strips with small 6 (i.e., short, widely separated 
boundaries), respectively. In Eqs. (|A7p and (jA8p one has 



^-ttS/4 



T,pq{K) 



Y, L[(24'+4p-39)^-l]/12 _ (^ ^ _^) j (A9) 



and 



1= 



a = e-2-/^ 5,,^5,,((t)=^(t('--i)/24x 

r=2 

TTp . nq 



. TTpr , nqr 1 , . „^ . „y , .,_,_. 

X sm sm / sm — sm — , (AlO) 

3 4 / V 3 \)'^' 

respectively. The three functions I]pg(K) in Eq. (|A7p can 
be written as 



S,,(..) = «:[(4p-3,f-l]/12^+^„^^ 



n=l 



(All) 



where the coefficients a„ take the integer values 0, ±1 
as determined by ([A9p . For small a the three functions 
Spq{(7) have the following explicit forms: 



Sn = 21/2(t1/« + (T - (T^ - 21/2(74+1/8 _ ^5 



521 
522 



-21/ V/« + (7 - (7^ + 21/^4+1/ 

„ „2 I „5 I 



(A12) 
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For closely spaced boundaries with i ^ j the deriva- 
tive with respect to k in Eq. (|A7p contributes due to the 
prefactor in Eq. (jAlip even in leading order W/L ^ 1 so 
that A.ij(O) depends on [i,j) as given below Eq. (|2.9p . 
For widely spaced boundaries (i, j) the derivative with re- 
spect to a in Eq. (|A8|) does not contribute to the leading 
behavior and 

A, j{l/5) -^ 7r(5"Vl2 = -Ap(0)/5^ 1/J ^ oo , 

(A13) 
is independent of («,j). In this case of W/L ^ 1 the 
stress tensor averages of the strip ST arc dominated by 
the periodic boundary condition in u direction and deter- 
mined by (T,, ||)st = (r„„)sT ^ Ap(0)/L2 with Ap(0) 
given below Eq. (|2.18p . 

The above expressions in Eqs. (|A5|) - (|A12p do not only 
serve to provide the aspect ratio dependence of the lead- 
ing contribution lS.Fi in Eq. (|2.17p of the free energy re- 
quired to rotate the small needle in the strip, but also 
to calculate the effective interactions for certain parti- 
cles of arbitrary size: in the first entry of Ref. [2^ for 
two particles of circular shape as well as in Appendix IB 21 
below for certain configurations of two needles in the un- 
bounded plane, of one needle in the half plane, and of 
one needle in a strip. 

Now we present the explicit forms of the scaling func- 
tions fQ'^\vN/W,{)) of the density profiles in Eq. (PTT)) 
in a strip of infinite length i = oo. These can be inferred 
from, e.g., Ref. |49j . In our notation they are given by 



ir^' = , /i°'°) = C/2 

4+'+) = (2C)V8 , /(+■+) = -C/2 



/(+•-) = -(2C)i/S , /(+'-) = (C/2)[3 - 4^2] 



/(+'°) = (C/2)i/8[i-,]i/4 



, /i+'O) = (C/2),s 



where 



C = 71/ cos{t:vn/W), s = sm{TTVN/W) . (A15) 



One can easily check that near the boundaries the cor- 
responding profiles {0{r]\f)) reduce to the half plane 
limits (©(rjv)) half plane with the amplitudes provided in 
Eq. (|2.5p . In particular, in the strip with i = + the three 

profiles f[ ' given in Eq. (jA14p exhibit the behavior 

fi'''\vN/W ^ -1/2) ^ 2^/^iivN/W) + (l/2))-V8 = 

{2W/aNy^^ corresponding to the half plane with bound- 
ary +. 



Appendix B: NEEDLES OF ARBITRARY 
LENGTH 



1. Symmetry preserving cases 

Needle and both boundaries belonging to O 



Here we establish Eqs. p.2p and (|3.13p for the free 
energy associated with the insertion of a needle in the 
case that both the needle and the boundaries are of the 
symmetry preserving "ordinary" type. We start with 
Burkhardt's result [69| for the thermal average of the 
stress tensor induced by n nonoverlapping "ordinary" 
needles embedded in the x axis of the unbounded (x, y) 
plane. If the n needles extend from a;i< to a;i>, from 
X2< to X2>, ... , and from a::„< to a:„>, respectively, with 
arbitrary real numbers a;i< < a;i> < X2< < a;2> ^ ■•■ < 
Xn< < Xn>, the stress tensor averages {Tki{x,y)) at a 
point (a;, y) follow from the analytic function 



{T{z))^io][o]...[o])^2-6f L 



1 



Z - ^1<; Z — 37^^ 



Z - X2< Z - X2> 



■2 •^n< ^ ■^n> 



(Bl) 



of the complex variable z = x + iy via [7l| the rela- 
tions {T^^{x,y)) = -{Tyy{x,y)) = -Rc(T(z))/7r and 
{Txy{x,y)) = {Tyx{x,y)) = lia{T{z))/'K where Re and 
Im denote real and imaginary parts, respectively. We 
point out the consistency of Eq. (jBip in the special cases 
Xm< -^ Xm>, hi which the needle m disappears, and for 
Xm> -^ Xm+K, in which two consecutive needles m and 
m+l merge into a single one. For n = 1 this reproduces 
the expression {T{z))^^'^^^ which is independent [40] of h 
and follows from mapping the half plane onto the entire 
plane outside a single needle (cp. Eqs. (A8) and (A9) in 
Ref. m. 

In order to derive Eq. p.2p we use Eq. (jBip for n = 2, 
put (a;i<,a;i>) = (— cx),0), so that the first needle fills the 
negative x-axis, and denote (a;2<,a;2>) by (a;<,a;>). The 
f A14) conformal transformation z = c^ /I with c = a + ib maps 
the {x, y) plane with the two "ordinary" needles onto the 
half plane (a > 0,6) with the "ordinary" boundary line 
a = and a single embedded "ordinary" needle extending 
from (a = y/lx^ =a<,6 = 0)to(a = y/Tx^ = ay,b = 0) . 
This leads to a stress tensor function 



where 



(Tic)) =~A/8 + c^A^/l6 



A 



(B2) 



(B3) 



In order to obtain this result one uses the transforma- 
tion formula for the stress tensor which includes the 
Schwartz derivative |7l| . The arbitrary length /, intro- 
duced for dimensional reasons, does not appear in the 
relation between {T{c)) and a<:,ay. Shifting the nee- 
dle away from the boundary line, i.e., increasing the dis- 
tance ajv = (a> +a<)/2 of its midpoint from the bound- 
ary while keeping its length a> — a< and orientation 
fixed, leads to a change in the free energy F± determined 
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by 



d F_L 



db{T^^{a,h)) 



(B4) 



where {T^^{a,b)) = — Re(r(c))/7r is the diagonal stress 
tensor component perpendicular to the boundary line 
a = 0. The integration path must extend between the 
boundary and the needle, i.e., < a < a< in Eq. (JB4|. 
In this region the integral over b is independent of a 
and is carried out most easily for a = 0. Integrating 
the result with respect to ajv and denoting the needle 
length a> — a< by Z) leads to the result in Eq. p.2p for 



In order to establish Eq. p.l3p one uses the transfor- 
mation c/l — exp(7rw/VF), w = u -\- iv, iw order to map 
the half plane (a > 0, b) with its embedded needle onto 
a needle in strip geometry as described in the context of 
Eq. (|37T3)) . With the ends of the needle at -u = ±D/2, 
this way one finds the stress tensor function 



{T{w)) = (7r/M^)2(-l/48 - S/8 + B"^ /l&) (B5) 



where 



B = 



sinli 9 



coshfl — cosh(27rii;/W^) 



(B6) 



with 6 defined in Eq. p. lip . The free energy change 
{d/dW)F^. ' ' /{ksT) upon widening the strip follows 
from the right hand side of Eq. (jB4p by replacing (a, b) 
by (u, w) and identifying (T1j.(m, w)) with Re(T(w))/7r. 
Here and below Eq. (|B4p the real parts of {T(w)) and 
{T{c)) enter with a plus and minus sign, respectively, be- 
cause in the complex w and c planes the directions _L per- 
pendicular to the boundaries point along the imaginary 
and the real axis, respectively. Performing the integral 
and integrating with respect to W leads to the result for 
F||°[°l°^ given in Eq. (1X^1) . 

Needle in a strip with periodic boundary condition 
Now we consider the geometry corresponding to Eq. 
p.l6p . For the infinitely long strip with periodic bound- 
ary condition containing a needle in parallel direction the 
stress tensor function {T(w)) follows from Eq. (|B1P with 
n = 1 and the transformation z/l = exp(27rz«/VF) yield- 
ing 



{T{w)) = (VW^)2(-l/12 + SVl6) 



(B7) 



with B given in Eq. (|B6p . Proceeding as for the [0,0) 
strip above, Eq. (|B7p leads to the result for F|| given in 
Eq. dSUl). 

Ordinary needle and + boundaries 

Finally, we derive Eqs. p.3p and p.l4p from the ex- 
pression |69[ 

(r(z))([^i[oi) = 2-6x 

I 1 1 I \ 



Xl< 



Xly 



X2< 



X2> 



for the stress tensor function induced by an "extraordi- 
nary" needle E extending from xk to xiy and an "or- 
dinary" needle from X2< to X2> ■ An "extraordinary" 
needle preserves the (-1- O — ) symmetry and can be 
realized in a lattice model as a line of spins with in- 
finitely strong nearest neighbor ferromagnetic couplings 
between them so that they all point either in the -I- or 
all in the — direction [3^ |4J]. The partition function 
Zmm = zil+m) + ^d-Hoi) in the presence of the 
two needles E and O differs by only a factor of 2 from 
the two identical partition functions Z^I+H^D = Z([~I[°1) 
in the presence of two needles + and O or — and O. 
(This argument holds if the partition functions are finite. 
This can be achieved by enclosing the whole system in 
a large box with "ordinary" boundaries.) Thus the cor- 
responding free energies differ by an additive constant 
—ksT ln2 which drops out from the free energy differ- 
ence upon changing the needle geometry as well as from 
the stress tensor so that 



{T{z)) 



{[EllO]) _ 



{T{z)) 



{[+][o]) _ 



{T{z)) 



{[-l[o]) 



(B9) 



Note the different sign sequences in (r(z))([-^l['^l' on the 
rhs of Eq. dnil) and in (r(z))([01[0]) ^j^ tj^g j-j^g ^f ^q. (|BT1) 
with n = 2. This implies that for obtaining the (+[0]) 
case Eqs. (|B2p and p.2p have to be modified by replac- 
ing ^ — > — A and ?? — > —i), respectively, which leads 
to E q dSin D- Similarly, for the {+[0]+) case Eqs. ([BS)) 
and p.l3p have to be modified by replacing B -^ —B 
and 6 — >■ —6, respectively, in order to obtain Eq. p.l4p . 



2. Cases w^ith genuinely broken symmetry 

In this subsection we derive the general expressions in 
Eqs. p.4p and (|3.15p for the Casimir forces which encom- 
pass also the cases {+[+]), (-[+]) and (+[+]+), (-[+]-) 
which cannot be reduced to cases with symmetry preserv- 
ing boundaries. We start by considering two needles 1 
and 2 on the x axis, similar as in Appendix IB 11 but with 
arbitrary universality classes [i] and [j] , respectively. The 
corresponding stress tensor (T(z))([*l[-'l'' can be inferred 
from the difference 

(T(z))(Hb1)_(r(z))([01[01) = 
1 



(z - Xi<)(z - Xi:^){z - X2<){z - X2>) 
{Xly -XI^){X2> -X2<) 



{i-kY 



where 



nAk) - -(-l + 6fc-fc2) 



-Ti J (fc) 



{2K{k)f 



(BIO) 



(Bll) 



(B8) 



and (r(z))'^t'^l['^l' is the stress tensor for two "ordinary" 
needles given by Eq. (|Bip with n ~ 2. The quantity 
< A: < 1 is related to cross ratios of the four needle end 
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points via one of the two following equivalent relations: 

(Xl> - X1<)(X2> - X2<) (1 - fc)^ 



and 



(x2< - xiy){x2> -a;i<) 



4fc 



and 



{x2< -Xi<_){x2> -a:i>) _ (1 + fc)^ 
(a;2< -a;i>)(x2> -xk) 



4fc 



(B12) 



(B13) 



The amplitude functions Aij{l/5) are given by 
Eqs. ([M]) - (|A10|) with the argument 1/(5 related to k 
via 



1/5 = A'(fc)/A'*(fc) 



(B14) 



with the complete elliptic integral [52[ functions K{k) 



and K*{k) = K{Vl-k^) (see Eq. (IXSll V 

Note that the pole of second order at z equal to a 
needle end, present with equal residues 2~^ in the stress 
tensor averages in Eqs. (|B1I) and (|B8p . is absent in the 
difference of averages in Eq. (jBlOp . This implies that, 
e.g., near the needle end at a;i< the leading contribu- 
tion to the corresponding average of the stress tensor el- 
ements Tki (see below Eq. (|B1[) ) is independent not only 
[4q | of the length and universality class of needle 1 but 
also of the presence (i.e., distance, length, and universal- 
ity class) of another needle 2. Moreover - as shown by 
Eqs. ((B2)) (|B5]) . and ([BZl above as weh as by Eqs. (|BT7)) 
and (|B28p below - this contribution is also independent 
of the presence of the concomitant boundaries of the half 
space and the strip. 

Equations (jBlOp - (|B14p follow from the Schwarz- 
Christoffel transformation |6j, UM which conformally 
maps the z = x + iy plane with the two needles [i] and 
[j] embedded in the x axis onto the rectangle or strip ST 
with boundaries {i,j) at -u = ±W/2, periodic boundary 
condition in u direction, and Casimir amplitudes A^.o, as 



introduced in the paragraph containing Eqs. (j2.7p - (|2.9p . 
The Schwartz derivative in the corresponding transfor- 
mation law of stress tensors [33, [71[ drops out from the 
stress tensor difference given by Eq. (|B10[) . 

The vanishing of r^j in Eq. (jBlip for {i,j) = (0,0) 
provides, together with Eq. (JB14I) . another expression 
for Ao,o{^/S) besides (but equivalent to) the ones in 
Eqs. dMI, dSl), and (jM]). For A+,o(l/'5) the cor- 
responding other expression follows from Eqs. (|B10[) 
and (|B11[) with {i,j) ~ {+,0) when combined with 
Eqs. dHIl), jHSl, and (lB9|) . 

The difference of the stress tensors in Eq. (jBlOp must 
vanish in the limit of distant needles with no correlation 
between them because the stress tensor for a single nee- 
dle in unbounded space is independent of its universality 
class [401. -f'o^ the same reason it must also vanish if 
two needles of the same universality class i = j ~ + (or 
— ) come close so that they merge, i.e., a;i> = X2<, and 
form a single -I- (or — ) needle. These expectations are in 
agreement with the behaviors 



liuik^i Tij{k) = 



(B15) 



limfe^oT,,^(fc) = [A»,»(0) - A,,j(0)]/7r 



(B16) 



of Tij for distant and close needles with S \ 0, k /^ 1 
and S /^ oo, k \0, respectively. Equation (|B15[) follows 
from Eqs. (|XT3)) and ([BU)) . 

Proceeding as described above Eq. (|B2p . the difference 
of the stress tensors for the needle in the half plane is 
obtained as 



(Tic)) 



(■'[j]) 



{Tie)) 



mo]] 



AA 



(1-fc) 



nAk) (B17) 



2 'J 



with (T(c))(°[°l) = (r(c)) from Eq. ([B2|) and the func- 
tion A from Eq. ([B3|. Via 



a</ay ^2k^^^/{l + k), 

fe=(l-l?l/2)V(l+l?l/')' 



(B18) 



k is related to the needle parameters a^/a^ and i? — 
(a> — a^)/{ay -f- a<) introduced above Eq. (jB2p and 
in Eq. (|2.2ip . Using the relation between k and i? in 
Eq. (|B18p together with the suitable functional relations 



K{k)^{l+d^l^fK*{^)/A, 
K*{k) = {l + d^''^fK{'d) 



(B19) 



between ellip tic integrals (see Eqs. 8.126.1 and 8.126.3 
in Ref. [S^]) yields together with Eq. ([Blip the more 
convenient expression 



{i-kY 



2'^^3ik) = ^nj(^) 



(B20) 



for the amplitude in Eq. (jB17p where f^ j is taken from 
Eqs. ([3?5 |) - ([3?7| . Equation ([BJ]) finally yields the expres- 
sion 



g ^(^bl) _ piO[o]) 



daN 



keT 



(B21) 



for the difference of the Casimir forces acting on the nee- 
dle in the cases (i, j) and (O, O) which implies Eq. p.4p 
upon renaming the dummy index j as h. 

For the special combinations (4-,0), (O, +), (— ,0), 
and (O, — ) the expression in Eq. (|B20p becomes indepen- 
dent of k and d and equals —1/4. This follows from com- 
paring Eq. (jBlOp with the simple expressions in Eqs. (|Bip 
and (jB8p and by using Eq. (|B9p . For combinations (i, j) 
with both i and j being "normal" universality classes the 
expression in Eq. (|B20p displays a nontrivial dependence 
on k and d. 

In Eq. p.6p the dependence on d of Ti_h and pi.^,, with 
(i, h) being arbitrary, can be readily calculated for a nee- 
dle nearly touching the boundary (J) = \/l — i?^ \^ 
0, (5 /^ oo) as well as for a distant or small needle 
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(^ \ 0, (5 \ 0) by combining the first and second ex- 
pression in Eq. ((XB)) with Eqs. (|T7)) and (jAlip and with 
Eqs. (jA8l) . ()Al2j) . and ()Al3p . respectively. In obtaining 
these limiting behaviors one uses the fact that Eqs. (|A9[) 
and (|A10|) and Eqs. ((3J)) and ((3^ imply the simple de- 
pendences 



K. = [(?9/4)P(r ) 



2m2 



(?9/4)P(792) (B22) 



of K and (T in Eqs. (|A7|) and ()A8|) on i? and ^, where 
P{x'^) can be expressed in terms of a power series in x^ 
given by 



P{x^) = exp{-Q(x^)/i?(a;2)} = 1 + x^/A + 0{x^) 



with 



(B23) 



Rix^) = {2/tt)K{x) = 1 + x'^/A + Oix'^) (B24) 



and (see Ref. [55 



Q{x^) = K{^/l-x'^) - R{x^) ln(4/a;) 

2 

1 , 




1x2 



1 X 3 
2x4 



1 



1 



1x2 3x4 



(B25) 



We use these relations to show that for small d the 
exact expressions in Eq. p.4p for the Casimir forces are 
consistent with the "small needle expansion" . To this end 
we first note that expanding the expressions in Eq. p.4p 
yields 



d F 



i^VA) 



ajv 



daM ksT 



-'i'^^^ 



0(1?") (B26) 



with the curly bracket being identical to the curly bracket 
in the expression for A^ ;, in Eq. (jASp . This follows from 
inserting Eqs. p.6p and (|A8[) into the second expression 
in Eq. (|3.4p and by using Eq. (|B24p for K{i}) as well as 
the relati on ad/ da = (1 - ^^ /2 + 0{'d^))dd/dd due to 
Eq. (jB22p . Furthermore, apart from terms of order d'^ the 
curly bracket in Eq. (|B26[) is equal to 1 + C/ + Ca with the 
expressions Q and C,a from Eq. ()2.19p for a small needle 
in a half plane. For example, in the case {i[h]) = (—[+]), 
with the small needle approximation |4l| 



-flAT 



d F 



(-[+]) 



dttN ksT 
d 



-^— ln(l-fC/ + a) 



-^— Infl - 2i/'*i9^/*^ + ?9/4 
dv 

-2i/4(3/32)i92+i/8 + (5/43)^3^ ^ 

(B27) 

the aforementioned relation between the curly bracket 
and l + (j + (a follows from the second part of Eq. (|A12[) 



and from Eqs. (jB22|) and (|B23p . Accordingly, the differ- 
ence between the exact expression for the critical Casimir 
forces in Eq. p.4p and their small needle approximations 
is of the order "j?^. 

Proceeding as in the paragraph containing Eqs. (jB5P 
and (|B6p , for the geometry of a needle of class j extend- 
ing along the midline of an (i, i) strip, as discussed in 
paragraph (ii) in Sec. IIIIl one obtains 



(r(y;))(«b]») _ (TH)(o[oio) = -Bx {tt/WY' 



i..(i)/(4i) 

(B28) 

{0[0]0) ^ 



with the function B given by Eq. (|B6p . {T{w)) 
{T{w)) from Eq. ([BSl) . nj as in Eq. ((X5|) . and t related 
to 6* = ttD/W as stated in Eq. (jXTSj) . The force in 
Eq. p.l5p follows from the stress tensor difference in 
Eq. (|B28p upon replacing [j] by [h] and by using a so- 
called shift equation as in Eq. (jB4p . 

Similar as in the paragraph containing Eq. (|B26p . the 
expanded expression 

piilh]i) 



+c9^+^^s + de"^} + O{0^) . 



-(l/3)2-^ (B29) 



of the exact result in Eq. p.l5p with coefficients a — e 
follows from Eqs. dMl), dH]), (HH]), (iB22|) . and (|B23)) 
and is related to the "small needle approximation" (4l| 



-W 



d F. 



[tlh] 



dW ksT 

-9^7, ln(l + dO^'^ 
do 



^^in(i + o-a) 



fo 



2+1/ 



(B30) 



where C/ and ^^ are given by Eqs. (|2.13p and (j2.14p 
with needle center vn = at the midline of the (i, i) 
strip. Note that the term eO^ with e given in Eq. (|B29p 
equals the contribution {-k / 2){D / 2WY /^i,i with the crit- 
ical Casimir amplitudes A^^, = — 7r/48 for {i,i) strips 
which enter into Ca in Eq. (pT4)) and Eq. (|B30p . Unlike 
e, the other coefficients a, b, c, d, /, and g depend on the 
boundary and needle universality classes {i,h). The devi- 
ation of the exact expression —W[d/dW)Fl-,''^ ^''' /{ksT) 
from its "small needle approximation" is given by the 
difference of Eq. (|B29p and Eq. (|B30p and is of the order 
6*^ because 



ae 



f^be 



0. 



(B31) 



For example, in the case {i[h]i) = ( — [+]—), in which the 
expression for the argument of the logarithm in Eq. (|B30p 
equals Zii as given by Eq. (|5.4p . one has 



(B32) 



= -2i/«, 6=1/8, c = -(5/3)2-7+i/«, 
= -(l/3)2-^/ = -(l/3)2-5+l/^g = 
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which satisfy Eq. (|B3T|) . 

For the convenience of the reader we provide the ex- 
phcit expressions 



(/o[o]W)_ = -ln{l + W4)[l + (5/16)i?2]}, 
(/+[+lW) = -ln{l + 2i/'*^i/«[l + (3/32)^2] 

\ / sna 

(f-l+]W) = -ln{l-2i/^^i/«[l + (3/32)^2] 

\ / sna 



and 



{foio]o(^) 



hi 



i|l-(l/3)2-^6|2 

+ (0/8)[l-2-^02]|^ 

(/",,,_(.))_. -hr{l- (1/3)2- V 
-(20)i/«[l + (1/3)2-5^2] ^ (^/8)[^ _ 2-^02]| 

(B34) 

for the [41[ "smah needle approximations" (/,-j'^|)sna 

and (/l|,j^)_ of /4, ^ FJ^['^lV(fcBT) and f^,^^ ee 

F|[''''l'V('^-Br). For the cases (+[0]) and (+[0]+) the 
expressions foUow from those for (0[0]) and {0[0]0) by 
the replacements d ^ —d and 9^—9, respectively. 

Appendix C: GLOSSARY 

In order to ease the reading of the text, here we com- 
pile the symbols and notations used and explain their 
meanings. 



A : Expression given by Eq. (jB3[) which serves to 
display the dependence on c of the analytic functions 
(T(c))('[''l) (see Eqs. ^^ and (|BT7| as well as the text 
below Eq. (jB9p ) and thus to display the position depen- 
dence of the stress tensor averages for the geometries of a 
needle of boundary universality class h in the half plane 
with orientation perpendicular to the boundary of bound- 
ary universality class i. 

Aq : Universal amplitudes of the profiles 
(©(r)) half piano in the half plane with boundary univer- 
sality class h where O = (p and O = e are the normalized 
operators (see Eq. (j2.4p ') of the order parameter density 
and of the deviation of the energy density from its bulk 
value, respectively (see Eq. (|2.5p '). 



a and h : Half plane coordinates perpendicular and 
parallel to the boundary a ~{) (see above Eqs. (|2.5p and 

dia)). 

Oat : Distance of the needle center from the boundary 
of the half plane (or from the lower boundary of the strip) 
(see Fig. [2]). 

a< and a> : Distance of the closer and farther needle 
end from the boundary for a needle in the half plane 
with perpendicular orientation (see Fig.[5][b) and the text 
above Eqs. ((XT|) and ((B2]) ): ay ^ a<^+ D. 



B : Expression given by Eq. (|B6|) which serves to 
display the dependence on w of the analytic functions 
(r(u;)) (*[''!*) (see Eqs. (|B5l) and (|B28| as weU as the text 
below Eq. (|B9|)) and thus to display the position depen- 
dence of the stress tensor averages for the geometries of a 
needle of class h embedded in the midline of a strip with 
(i, i) boundaries. 

c = a + ib : Complex variable specifying the position 
vector {a,b) in the half plane (sec above Eq. (|B2[) '). 

Qj = [AFni/{kBT)]/{D/W)'^ : Normalized next-to- 
leading contribution to the quasi-torque acting on a small 
"ordinary" needle O in an (i,j) strip (see the paragraph 
following Eqs. (|5.2p and (|5.3p and the caption of Fig. [TO|) . 



D : The number of missing bonds (fixed spins) in the 
lattice description of the "ordinary" ( "normal" ) needle in 
Sees. UandUVKsee Fig. [Das weU as Figs. Hand El) or ^ 
the length of the needle in the continuum descriptions 
used in Sees. HI HIl |llTl and Appendix [B] (see Fig. (2). 

E : "Extraordinary" boundary with infinitely strong 
ferromagnetic nearest neighbor couplings between sur- 
face spins (see Ref. [38j). 

i^ll and F± : Free energy cost to transfer the needle 
from the bulk into the strip (or into the half plane) with 
its orientation parallel and perpendicular to the bound- 
aries, respectively (see Eqs. (|2.10p - (|2.14p and the text 
above Eq. p.ll[) or below Eq. p.ip ). In particular, F^\ = 



F,p[''l') = kBTfl,^{9) and F^ 



F 



iz[h]) 



f^BTf^^^y 



for needles h of arbitrary length embedded in the mid- 
line of an (i,i) strip and perpendicular to the boundary 
i of the half plane, respectively. 

AF = F|| — F±_ = Fi — Fq : Free energy required to 
turn the needle about its center from the perpendicular 
to the parallel orientation (see Eq. (|2.15p V 

AFi and AF„; : Leading and next-to-leading contri- 
butions, respectively, to AF for the cases («[0]j) and 
{0[+/-]0) of a small needle (see Eqs. (jri^ - ipi^ 
as well as the corresponding half plane relations in 
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Eq. (12:201) ). 



K 



(«j) 



ST 



Free energy of the stri p S T without needle and 
with boundaries (i,j) (sec Ref. [43|). 

Fcr(A) : The free energy belonging to the lattice 
crossover Hamiltonian (see the paragraph containing 
Eqs. ga-lMl))- 



h : Characterizes the universality class of the needle 
surface. 

{i,j) : Characterizes the surface universality classes of 
the (lower, upper) boundary of the strip. 

{i[h]) : Needle of class h in the half plane with bound- 
ary of class i. 



M[h\i) 



fl\^^^{9) = Fp'''^'>/{kBT) : See Eq. ^A^ as wcU as 
Figs. SI and El 

f^f^i^) = Ff^''^^/{kBT) : See Fig. El and below 
Eq. dtH). 



/b and /s : Bulk free energy per area and surface free 
energy per length (see Ref. |42|). 

Pq : Universal scaling functions of the profiles 
(C'(r))sT for O = 0, e in the needle-free strip ST with 
boundaries (i, j) (see Eqs. (|2.7p and (|A14p ). 



{i[h]j) : Needle of class h in the strip with boundaries 
of classes i and j . 

i/ > : Ferromagnetic coupling strength between 
nearest neighbor Ising spins on the square lattice im- 
plying the bulk critical temperature Tc = [2/ln(-\/2 -I- 
l)](J'/fcB) (see Eq. (|4.ip and the paragraph above Eq. 

63)). 



Ju,v;u' ,v'J ■ Coupling between nearest neighbors {u, v) 
and (m'j v') in a strip containing a needle of broken bonds 
(see Fig.ElJa) andjHl^b)). J equals and 1 for broken and 
unbroken bonds, respectively. 



(p) 
/e : Position-independent universal scaling function 

of (e)sT in the needle-free strip ST with double peri- 
odic boundary conditions (see below Eq. ()2.18|) and Ap- 
pendix |XT]). 



k : Parameter characterizing the configuration of two 
needles embedded in the x axis of the unbounded (x, y) 
plane via the cross ratio of their endpoints (see Eqs. (|B12[) 
and (IBT311 '). 



9i\h] (^) ^^'^ 5!Lii(^) • Scaling functions for the effective 
force acting on a needle h perpendicular to the boundary 
i of the half plane and for the disjoining force induced in 
an (i, i) strip by a needle h embedded in its midline (see 
the text below Eq. (|3T|) or above Eq. (|3ll|) ). 

"HsT : Lattice Hamiltonian for strips without a needle 
(see Eq. gJ])). 

•HsT + H^]^^ = Ha and Hst + uf^ = Ui : Lattice 
Hamiltonians for strips containing an embedded needle 
with orientation perpendicular and parallel to the bound- 
aries, respectively (see Eqs. (|4.2p and 

A 17 — 17 17 — njU^) 

Hamiltonians for parallel and perpendicular orientation 
of the needle h with a fixed center (see Eqs. (|4.4p - (|4.7l) '). 



'W^ ' : Difference of lattice 



/(O) 



'Hc!^''(A) and Her (A) : Crossover Hamiltonian for the 
needle of broken bonds and of fixed spins, respectively 
(see Eqs. (|iTUl) and (pi^ as weU as Figs.i^c) andlUJc)). 



H'^^^ and 'H}-^^ : Lattice Hamiltonian for a strip con- 
taining a cross-shaped hole with the two bars correspond- 
ing to the two orientations of the needle of broken bonds 
and of fixed spins, respectively (see below Eqs. (|4.10p 
and gm). 



L : The number of columns in the lattice model for 
the strips considered in Sees. |I] and IIVI (see Fig. [T] and 
Figs.|6l[8]) or [23| the length of the strip in the continuum 
descriptions used in Sccs.HIIh IIHH and Appendix lAl (see 
Fig. El). 

I : Arbitrary length in the conformal transformations 
given in the text above Eqs. (|B2|, HI]), and dBT]) which 
relate different geometries. It is introduced for dimen- 
sional reasons only and, due to dilatation invariance, 
drops out from equations relating quantities which be- 
long to the same geometry. 

n : Unit vector describing the orientation of the needle 
(see Eq. ((O)) ). In the strip, n = (?i|| , n^) (sec Eq. ([^ l. 



O : The "ordinary" surface universality class. Cor- 
responding boundaries or needles induce disorder in the 
system of Ising spins, i.e., in their vicinity the probability 
to find parallel nearest neighbor spins is smaller than in 
the bulk. 

-|- and — : The two "normal" surface universality 
classes with the tendency to order the Ising spins in the 
+ and — directions, respectively. 

rpf : Position vector of the center of the needle (see 
Eqs. (|2.2p and (|2.3I) '). In the strip, r^ = {un,vn) (see 
Eq. JSIM ). 
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Si and 5*^ : Operator contribution to the normalized 
Boltzmann weight of the small needle which is isotropic 
and anisotropic, respectively, with respect to the orien- 
tation of the needle (see Eqs. (PTT|) -(P3 |) '). 

ST : Denotes the strip in the absence of the needle. 



Z'li and Z± : Partition functions corresponding to F|| 
and F± (see Eqs. H^JU^ and (P1T|) ). 






Partition function of a M^ x L strip ST without 



needle and with boundaries (i, j) (sec Refs. [43 and \44 



"sna" ("small needle approximation"): Truncated 
form of the "small needle expansion" (see the paragraph 
containing Eqs. (|2.16|) - (|2.18|) and the first paragraph of 
Subsec. IV Cp . For explicit expressions see Eqs. (|5.4p . 
(IB331) . and (lB34)) . 

So : Exterior spin fixed to the value +1 (see the first 
paragraph in Sec lIVB^ as well as Fig. Hl^c)). 

Tki : Stress-tensor operator (see Eq. p.Sp ) with its el- 
ements in the strip denoted by Tu 
(see Eq. (1^ ). 



Til 



r^||, and Tj_j_ 



{T{z)), {T{c)), and (T(w)) : Analytic functions in the 
unbounded plane with needles, in the half plane with 
a needle, and in the strip with a needle, respectively, 
which determine the corresponding stress tensor averages 
as explained in Appendix IB II (see Eqs. ([BT]) . (|B2|) . (|B5t . 
(iBTll- dBTOl) . (IBTTI) . and (lB28l) as weU as Ref. (Zll)- 



t = tanh(6'/2) : Useful short notation according to 
Eqs. (|XT5)) and ((B28)) . 

u and V : Strip coordinates j27| parallel and perpen- 
dicular, respectively, to the boundaries of the strip. In 
the continuum description the boundaries of the strip are 
at u = ±W/2 (see Fig. [5]). In the lattice description of a 
"normal" ( "ordinary" ) needle the coordinates u and v of 
the lattice vertices have integer (half odd integer) values 
(see Fig. [Hand Figs. mi). 

uat and vjy : Coordinates parallel and perpendic- 
ular, respectively, to the strip of the position vector 
rj\[ = {un,vn) of the center of the needle. As explained 



in Sec. [J and Ref. [27[, un and vn are lengths in the 



continuum description used in Sec. Ulland in Eqs. (|A14[) 
and (JAISP of Appendix [XI while in the lattice description 
used in Sec. IIVI thev arc measured in units of the lattice 
constant and have integer values for both "ordinary" and 
"normal" needles. 

W : The number of rows in the lattice model for the 
strip or [27(] the width of the strip in the continuum de- 
scription. 

w = u + iv : Complex variable specifying the position 
vector {u,v) in the strip (see above Eq. ([B5|). 

x^ = 1/8 and x^ ~ 1 : Scaling dimensions of the order 
parameter and energy densities, respectively (see below 
Eq. ([H])). 



^(['ii]['i2]) . Partition function of a large system con- 
taining two needles hi and /12 (see below Eq. (jB8[) ). 

z = X + iy : Complex variable specifying the po- 
sition vector {x,y) in the unbounded plane (see below 

Eq. m\i ). 



Aij(l/6) : Casimir amplitude describing the univer- 
sal contribution ~L^^d<^''Q^\S)/dW = A,^j{l/ S)/W^ to 
the disjoining pressure per fc^T of an {i,j) strip without 
needle (see Eqs. ([^ . ([^ and Appendix [X2|) . 



Aij = Aij(0) : Casimir amplitude for an (i, j) strip 
of infinite length L = 00 and without needle (see below 
Eq. l[Tm ). 

Ap{l/6) : Casimir amplitude for the double periodic 
strip without needle (see below Eq. (|2.18p and Appendix 
IXTIl . 

Ap = Ap(0) = -7r/12. 

S = L/W : Characterizes the shape of the strip (see 
Eq. ([M])). We caU W/L = 1/(5 the aspect ratio of the 
strip. 

e(r) : Energy-density operator with its average in the 
unbounded plane (bulk) at bulk criticality subtracted 
(see below Eq. (|2.3p i and normalized according to Eq. 



C/ and C^A '■ Contributions to Zn and Z± which arise 
via the "small needle expansion" from the operator 
contributions Sj and Sa to the normalized Boltzmann 
weight of the needle; they arc isotropic and anisotropic 
with respect to the needle orientation, respectively (see 
Eqs. ([23D-([23D, (irra . and SUM ). 

d = D/{2aN) : Characterizes the size versus the dis- 
tance to the boundary for a needle in the half plane with 
its orientation perpendicular to the boundary line (see 
Eq. ([S3D). 



d = Vl — ^^ = ^ci<a>/aN ■ Approaches zero if the 
closer end of the needle approaches the boundary. 

■i^o ■ Threshold value of 1} above which the interac- 
tion /^r,i(^) between a — boundary of the half plane 
and a -I- needle perpendicular to it deviates significantly 
from the corresponding "small needle approximation" , 
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i.e., from the last equation in Eq. (jB33|) (see the discus- 
sion of Fig. EJa) in the paragraph between Eqs. p.lOp 
and dSH])). 

6 = ttD/W : Characterizes the size versus the dis- 
tance to the boundaries for a needle embedded within 
the midline of the strip (see Eq. p.!!}) ). 

K = exp(— 7r(5/4) : Variable used for the aspect ratio 
dependence of A^ ^ in the case of a long strip with L/W = 
(5^1 and close boundaries (see Eqs. (|A7|) and (jA9p ). 

A^ and A^ : Strengths of the coupling to the lower 
and upper additional outside row, respectively, of fixed 
spins generating strip boundaries i and j of "ordinary" 
or "normal" character in the lattice model (see Eq. (|4.1|) 
and Figs, m and Umi). 

A : Parameter within the lattice model describing the 
crossover of the needle orientation from perpendicular 
(A =: 0) to parallel (A = 1) orientation with respect to the 
boundaries of the strip (see Eq. (|4.7p as well as Figs.[ni[c) 
andlHIc)). 

Pi^h{'&) '■ Auxiliary function determining the Casimir 
force on a needle h in the half plane with boundary i and 
determining the disjoining force induced in an (i, i) strip 
of infinite length upon inserting a needle h (see Eqs. p.4p 
and ((3TT5)) . respectively). Via Eqs. ^^ and ^J^ the 
function pi^h{'&) is related to the dependence Ai_ft,(l/(5) on 
the aspect ratio of the Casimir amplitude of the needle- 
free strip with boundaries (i, h). 

^ and ^ : Sum of products of those nearest 

(inc.) (deer.) 

neighbor spins in the crossover Hamiltonian "Her (A) the 
ferromagnetic coupling strength of which increases and 
decreases, respectively, upon increasing A (see Eqs. 
(PJU)) andFig.ElJc)). 



(+) 

^ : Sum of those four spins which are coupled to 

(one) 

an external spin (and carry a northeast arrow) for both 
needle orientations shown in Figs. [SJa) and [Hl[b) (see 
Eq. glU)). 

^ : Sum of products of the center spin of the "nor- 

(zcro) 

mal" (-I-) needle and of its four nearest neighbor spins 
(see Eq. (|4.13p ). These four products are missing in the 
crossover Hamiltonian "Her (see Eq. (|4.19p ) and the four 
corresponding couplings are absent in Figs. [Slja), [5Jb), 
andilc). 

^ and ^ : Sum of products of those near- 

(inc.) (deer.) 

est neighbor spins in the crossover Hamiltonian Her (A) 



the ferromagnetic coupling strength of which increases 
and decreases, respectively, upon increasing A (see 
Eqs. dHHl), dmni), and (|ilt| as well as Fig.HlJc)). 

{+) (+) 

^ and ^ : Sum of those spins which in the 

(inc.) (deer.) 

crossover Hamiltonian Her (A) are coupled with increas- 
ing and decreasing strength, respectively, to the external 
spin So = 1 (see Eqs. (|iT7l) - (|i:Tg|) and Fig.HJc)). 

a = exp(— 27r/(5) : Variable used for the aspect ra- 
tio dependence of A^j in the case of a short strip with 
L/W = 5 <^ 1 and distant boundaries (see Eqs. (jASp 
and (ETOl) '). 

Ti,h{'d) : Auxiliary function with the same character as 
(and simply related to) Pi^h{'&) (see Eq. p.Sp ). 

'^ST v*^ = L/W) : Universal shape dependent and 
scale free contribution to the free energy per ksT of 
L X W strips with boundaries (i, j) but without needle 

(see Eq. (HH) and Ref. [H]). *gT^((5) denotes the corre- 
sponding contribution for the double periodic strip. 

(f>{r) : Order-parameter-density operator, normalized 
according to Eq. 



Xpqi^) or Xpq((5/2) : Auxiliary functions for the 
Casimir amplitudes of W^ x i strips without a needle and 
with double periodic boundary conditions or with (i,j) 
boundaries (see Eqs. (jAip or (|A6p as well as Ref. |43|). 

(...) : Thermal average which may be specified by 
means of subscripts and superscripts such as 

(•■•)buik : for the unbounded plane without embedded 
particles (see Eq. (|2.4p '). 

(■•■)haif plane : for the half plane without embedded par- 
ticles (see above Eq. (|2.5p '). 

(...)gT and (...)g!^-' : for a strip without embedded par- 
ticles and with boundaries (i, j) (sec Eqs. (|2.7p and (|2.8p ). 



(J _^^(i[/i]) . £qj. ^jjg YiQ}f space with boundary i and an 
embedded needle h (see Eq. (|B17p V 

(...)^*[''l-') : for the strip with boundaries (i,j) and an 
embedded needle h (see Eq. (jB28p V 

^^^^^([hi][h2]..[/i„]) . Jqj. ^Yie unbounded plane with n em- 
bedded needles hi, /12, •■, hn (see Eqs. (|Bip and (jBSp - 
(iBTOll ). 

(...)cr : Thermal average based on the lattice crossover 
Hamiltonian (see Eqs. (|4.5p and (|4.6p ). 
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parameter profile induced by the boundaries {i,j) there 
appear terms oc (D/Wy'^ inside the logarithm deter- 
mining A_F which cannot be made sufficiently small in 
simulations with mesoscopic D (in order not to exceed a 
tractable lattice size W x L) and prevent a useful expan- 
sion (see the remarks above Eq. (|2.16p ). Actually, for the 
case (—[-I-] — ) the naive application of (ii) would predict 
the wrong orientation of the needle near the strip cen- 
ter (compare the discussion in the paragraph containing 

Eq. dsaiD- 

[62] For the case (+[+] — ) the simulation data and the an- 
alytic prediction indicate that AF becomes positive, in 
agreement with the observation (i) in the Summary, for 
values vn/W which are closer to the + boundary (i.e., 
which are more negative) than the ones shown in Fig. 1151 

[63] See, e.g., V. I. Smirnov, A Course of Higher Mathematics 
(Pergamon, Oxford 1964), Vol. 3. 

[64] Unlike in d = 2, in three spatial dimensions (d — 3) the 
scaling dimension 2 + x^ of Srjy^ 9rjv, </> is smaller than the 
scaling dimension d of the stress tensor. Thus for small 
"normal" particles with uniaxial shape-symmetry (such 
as a disk of diameter D or a needle of length D) and lo- 
cated between infinitely extended parallel plates - a dis- 
tance W apart and with broken symmetry -, the leading 
orientation dependent contribution within the small par- 
ticle expansion is proportional to D'^'^^'^ times the second 
derivative of the order parameter profile in the system 
without the particle. Moreover with x^{d = 3) « 1/2 
(about four times larger than in d = 2) it might be possi- 
ble within a simulation to keep (D/W^)^* and (D/ajv)^* 
significantly smaller than 1 and to single out the leading 
quasi-torque AFi proportional to D^^^*. 

[65] Needles in d = 2 and disks in d = 3 are "impenetrable" 
boundaries of finite extent D in the sense that spins close 
but on different sides, separated by the width W of the 
boundary with W <^ D, are coupled only via correlations 
around the ends of the boundary. (This also applies to the 
needles in the present lattice model with nearest neigh- 
bor couplings and where W is of the order of the lattice 
constant.) Thus the density profiles they induce at meso- 
scopic distances much larger than W are independent of 
W and depend only on D. For needles in d = 3, however, 
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its normalization is specified by its transformation law. 
For two geometries G and G related by the conformal 
mapping z = z{z), the corresponding analytic func- 
tions {T{z))g and {T{i))^ are related by {T{z))q = 
{dz/dzf{T{z))G + S/2A with the Schwartz derivative 
S = {z'" /z') - {■i/2){z" /z'f where the primes denote 
derivatives with respect to z. For example, for the map- 
ping z = z{c) = c? /I considered above Eq. (|B2|) this re- 
lation allows one to obtain {T{c))q = {T{c)} in Eq. (|B2)) 
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[72] For the simple needle configuration a;i< — 

— 1/fc, Xiy = —1, X2< = 1, X2> = 1/fc with 

< fc < 1 the mapping is provided by 
u + iv = i{W/2)[J^ dt/^{l - i2)(l - k2t'^)]/K{k). 
This maps the upper and lower halves of the z plane 
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Crossing the x axis outside of the interval where the 
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which amounts to the periodic boundary condition 
in u direction. The more general needle configuration 
corresponding to Eq. HB10|I is related to the simple one 
by a special conformal transformation which preserves 
the real axis. 



